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ABSTRACT 

Extensive photometric studies of the globular clusters located towards the center of the Milky 
Way have been historically neglected. The presence of patchy differential reddening in front of these 
clusters has proven to be a significant obstacle to their detailed study. We present here a well-defined 
and reasonably homogeneous photometric database for 25 of the brightest Galactic globular clusters 
located in the direction of the inner Galaxy. These data were obtained in the V, and / bands 
using the Magellan 6.5m telescope and the Hubble Space Telescope. A new technique is extensively 
used in this paper to map the differential reddening in the individual cluster fields, and to produce 
cleaner, dereddened color-magnitude diagrams for all the clusters in the database. Subsequent papers 
will detail the astrophysical analysis of the cluster populations, and the properties of the obscuring 
material along the clusters' lines of sight. 

Subject headings: globular clusters: general - Galaxy: bulge - Galaxy: evolution - Hertzsprung- Russell 
and C-M diagrams - stars: horizontal-branch - dust, extinction 



1. INTRODUCTION 

The globular cluster system of the Milky Way has long 
been used to learn about the evolution of the Galaxy. As 
the Galactic globular clusters (GGCs) constitute some of 
the oldest systems in the Milky Way, the study of the 
stellar populations of these fossils can give us important 
clues of the early stages of the Galaxy's formation. The 
tool most extensively used in this task has been the anal 



ysis of the color- magnitude diagrams (CMDs) (e.g., de 
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Angeli et al .| ( [20051 ); [ Marm-Franch et al .| ([20091)). how- 
ever, the presence of significant differential extinction in 
low-latitude fields, particularly near the Galactic center, 
greatly complicates traditional CMD analyses. As a re- 
sult, the study of many GGCs located towards the inner 
Galaxy has been historically neglected. 

Various recent studies have tried to overcome the diffi- 
culties associated with the study of inner GGCs to better 
exploit them as probes of the stellar populations near the 
Galactic Center. One obvious approach is to use near- 
infrared photometry to study these clusters (e.g. Davidge 
( |2000D ; iValent i et all ( |2007D ) to take advantage ot the 
smaller extinction m these bands. But to extract precise 
information about the ages of the GGCs, an accurate 
location of the main sequence turn-off (M STO) point is 
crucial in most methods of cluster dating ( [Stetson et al. 



1996[ ISarajedini et al ||1997[ [Gratton et ar]|2QQ3D . In- 

trared photometry sumciently deep to reach the main 
sequence (MS) with enough precision to accurately lo- 
cate the turn-off (TO) point is difficult to achieve, and 
only now are we starting to see the ffist results after 
the careful applica tion of new adaptive optics techniques 
on big telescopes (Moretti et al. |2QQ9 ). Deep photom- 
etry for non-highly reddened GGCs is much more eas- 
ily obtained in the optical, but at the cost of increased 
extinction. A number of techniques have been devel- 
oped to produce extinction maps for individual clusters 
as a means of dealing with this issue (e. g., Mel bourne^ 
Guhathakurta] ( |2Q04D ; |von Braun fc M ateo (2UUll)j. in 
general, the resolution and accuracy of optical extinc- 
tion maps are not adequate to fully eliminate the effects 
of differential reddening at a level of precision to produce 
deep optical CMDs of similar high quality now typical for 
high-latit ude clusters with l ittle or no differentia l extin c- 
tion (e.g., [Rosenberg et al.|([2QQQaD;|Piotto et al.|(|2QQ2D) 



In this paper we present a new optical photometric 
database consisting of a sample of 25 GGCs located to- 
wards the Galactic Center. We use a new dereddening 
technique (Alonso- Garcia et aL] (|2011l), from now on re- 
ferred to as Paper 1) to map the differential extinction 
along their fields, and produce new, cleaner, differentially 
dereddened CMDs of the clusters. In section 2, we give 
the criteria used to define our sample of GGCs. Section 3 
summarizes the steps that we followed to obtain precise 
astrometry and optical photometry of the stars in these 
clusters suitable for our analyses. In section 4 we apply 
the dereddening technique described in Paper I, and pro- 
vide an overview of the differentially dereddened CMDs 
of the clusters in our sample, along with extinction maps 
along their fields. We also describe the characteristic fea- 
tures of the environments in which they are located, and 
briefly summarize previous optical and infrared photo- 
metric studies where appropriate. Finally, in section 5, 
we provide a summary of our results. 

2. SELECTION OF OUR SAMPLE MEMBERS 

The cluster sample presented in this study consists of 
25 GGCs in the direction of the inner Galaxy, all located 
within 30 deg of the Galactic Center. Because the preci- 
sion of our dereddening technique (Paper I) requires good 
sampling and photometry down to a few magnitudes be- 
low the MSTO, and also depends on the density of stars 
and the spatial variations in the extinction, the clusters 
in our sample were chosen to also satisfy the following 
criteria: 

• Exhibit moderate mean extinction, implying that 
they may suffer from extinction variations. We 
explicitly restricted our sample to clusters with a 
mean reddening of E{B — V) > 0.07 mag. 

• Be sufficiently luminous to possess a well-defined 
MS. Our analysis requires a significant number of 
stars (at least a few hundreds) to calculate the ex- 
tinction in a region. We therefore chose clusters 
with a luminosity satisfying My < —6. 

• Be relatively nearby. With our dereddening tech- 
nique, the stars that provide most of the infor- 
mation about the differential extinction are those 



located in the CMD sequences most nearly or- 
thogonal to the reddening vector (subgiant branch 
(SGB) and upper MS). Since one of the goals of 
this project is to accurately calculate the relative 
ages of these clusters, we must also reach the MS 
with good photometric precision at the TO in or- 
der to carry out a reliable age/metallicity analysis. 
Hence, we chose clusters with an apparent distance 
modulus of (m — M)y < 16.6. 

• Be sufficiently extended so that we can define maps 
that cover a significant solid angle around the clus- 
ters. We therefore chose clusters with a tidal radius 
Tt > 7.5 arcmin. 

We were able to observe 25 of the 31 GGCs that ful- 
fill these requirem ents, according to the 2003 version of 
the 



Harris 



( 1996 ) catalog. Their positions are shown in 

Figure [T] along with the rest of the clusters located in 
the inner Galaxy. Their characteristics, according to the 
most recent (December 2010) versiorj^of the Harris cat- 
alog, are given in Table [l] 

3. OBSERVATIONS AND DATA REDUCTION 

We obtained optical photometric data for our sample 
of GGCs using the Magellan 6.5m Baade Telescope lo- 
cated at the Las Campanas Observatory (LCO) in Chile, 
and the Hubble Space Telescope {HST). In this section we 
provide a complete summary of these observations and 
explain in detail the steps followed for the reduction and 
calibration of the data. The final photometry has an ab- 
solute precision with respect to the calibrating stars of 
a ~ 0.02 — 0.03 magnitudes in most cases. The internal 
precision among the stars observed in each field is around 
a ~ 0.01 — 0.02 magnitudes, sufficient to achieve the pre- 
cise extinction maps that will allow us to produce clean 
CMDs. From these we will aim to derive the cluster pa- 
rameters with an accuracy similar to those derived from 
high-Galactic latitude clusters CMDs that suffer little or 
no differential extinction. 

3.1. Observations 

The GGCs in our sample were observed over four 
nights. May, 30th to June, 2nd 2005, with the LCO 6.5 
m Magellan Baade Telescope, using the Inamori Magel- 
lan Areal Camera and Spectrograph (IMACS) in imag- 
ing mode. We used the f/4 camera to image a field of 
view (FOV) of 15.46' x 15.46^ with a pixel size of 0.11". 
All fields were obse rved usin g standard Johnson- Cousins 
5,V, and / ffiters (Bessell 1979). Two sets of observa- 
tions with different exposure times (short and long) were 
taken in the B and V filters, and three (extra-short, 
short, and long) in /, for every cluster. Although the 
nights during our observing run were not all completely 
photometric, the seeing conditions were very good, with 
average values of ~ 0.6^' in V for the whole run. Ta- 
ble [2] lists the details of these ground-based photometry 
observations. 

^ According to this updated version of the catalog, two of the 
clusters in our sample, NGC 6287 and NGC 6355, do not fulfill 
anymore the apparent distance modulus requirement to belong to 
our sample. However, we have still kept these clusters in our stud- 
ied sample, since previous studies of these objects have been very 
scarce, and their published parameters are accordingly quite un- 
certain. 
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In order to obtain useful photometry of the inner re- 
gions of the more centrally crowded clusters, we sup- 
plemented our Magellan images with images taken with 
the Advanced Camera for Surveys (ACS) aboard HST. 
These data were obtained in by our group's Snapshot 
program 10573. The ACS has a FOV of 3.37' x 3.3^, 
with a pixel size of 0.05'^ Five clusters of our sample 
were observed using the fA'ibw{B^^^)^ /555k;(V555), and 
/814k; (/814) filters. Table |3] lists the details of these new 
HST observations. 

To better calibrate our photometry, we also used im- 
ages available through the HST data archive for all 
the clusters in our sample. The data taken from 
the HST data archive are comprised of /439k; (^439), 
/555k;(V555), /606k;(F606)7 and /814k;(/814) images ob- 
tained with the Wide Field Planetary Camera 2 (WFC2), 
and of /435k;(5435), /606k;(F606)7 and /814k;(/814) im- 
ages taken with the ACS. Table [4| lists the different HST 
programs that we use data from. 

3.2. Data reduction and photometry derivation 

For the ground-based data, the initial processing of 
the raw CCD images was done with the routines in the 
ccdred package of the Image Reduction and Analysis Fa- 
cility (IRAF). The images were first corrected for bias, 
then flatfielded using a combination of dome and twi- 
light flats. Afterward, since different frames for every 
cluster were taken with different exposure times in the 
5, V, and / filters (see Table [2|, and they had small spa- 
tial offsets between them, thelrames of every individual 
cluster were aligned, using the IRAF task imalign^ and 
average- combined for every exposure time and filter, with 
the IRAF task imcombine. 

Stellar photometry was carried out on the ground- 
based pro cessed images using a n updated version of 
DoPHOT dSchechter et al.||1993[ ). This version works 
on any platform and accepts images consisting of real 
data values. It also provides better aperture corrections 
than previous versions, allowing for variations in aper- 
ture corrections as a function of field position and stellar 
magnitude (see Appendix [A]) . 

We also carried out an astrometric analysis of the clus- 
ter fields in our sample. We derived coordinates (Right 
Ascension a and Declination S) for all stars detected in 
our photometric study by comparison with bright stars 
obtained in each field from the Two Micron All Sky Sur- 
vey (2MASS) catalog stars available through the Infrared 
Processing and Analysis Center (IPAC) website. In prac- 
tice, more than 100 stars were available as astrometric 
references within the fields of every chip of our CCD cam- 
era. A third order polynomial fit, done with the IRAF 
task mscpeak^ produced dispersions of a ~ 0.25'^ con- 
sistent with the catalog precision. Using the astrometric 
information of the images we could also calculate if there 
were variations in the pixel area coverage accross the 
images, which could lead to miscalculations in the mea- 
sured fluxes. But the pixel area changed by less than 1% 
across any of the chips of the camera, and therefore we 
did not need to apply any corrections to the measured 
photometry. 

To calibrate our ground-based data, we first trans- 
formed the photometry from each of the individual eight 
CCDs in the camera in IMACS to a common zero point, 
and then transformed this system to a standard system- 



in our case, the Johnson-Cousins system ( Bessell|[l979 ). 
We arbitrarily chose one of the CCDs, chip 2, as the 
one defining the instrumental system and brought the 
photometry from the other chips to this system (see Ap- 
pendix IbI). After this was done, we placed our photom- 
etry in the John s on- Co usins photometric system by ob- 
serving [Landolt ( 1992 ) fields over a range of airmasses 
during the niglits of our 2005 observing run (see Ap- 
pendix [C|. Since these nights we re not completely pho- 
tometric we used Stetsonl (2000*) photometric standard 
stars and the photometryTrom the clusters obtained with 
the HST to fine-tune the calibration of our ground-based 
data (see Appendix [C]) . 

The photometry 01 the data from the HST observa- 
tions was obtained using the programs HSTPHOT for 
the WFPC2 data (Dolphin 2000) and DolPHOT for the 
ACS data ( jDolphin||2Q00p . These programs have been 
specifically tailored to analyze data from these cam- 
eras. They take the images for every cluster retrieved 
from the the HST in the different filters, and analyze 
them all at the same time. The programs align the 
images, combine them and provide the final combined- 
frame and individual-frame photometry, both in the HST 
filter system and in the Johnson-Cousins system. We fol- 
low the author's prescriptions to give the input param- 
eters for the analysis to both programs. Photometric 
goodness-of-fit parameters were employed to select only 
objects with high-quality measurements. Only measure- 
ments for which an object was classified as stellar (HST- 
PHOT and DolPHOT types 1 and 2), and for which 
—0.1 < sharpness < 0.1, crowding < 0.5, and errors 
in all measured magnitudes <j < 0.1 were retained. 

4. DEREDDENED COLOR-MAGNITUDE DIAGRAMS AND 
EXTINCTION MAPS 

In this section we present the results of applying the 
dereddening method described in Paper I to all the clus- 
ters in our sample of inner GGCs. We provide the CMDs 
of these clusters before and after correcting for differen- 
tial extinction and describe any improvements. We also 
provide the extinction maps that the dereddening tech- 
nique generates, and compare them w ith the extinction 
maps supplied by Schlegel et al. (1998), from now on re- 
ferred as the SFD extmction maps~~For each cluster, we 
provide references to its more recent - but not necessarily 
all- optical and infrared photometric studies. We also re- 
port studies that explore the presence of different stellar 
populations in the observed GGCs, and speculate on the 
possibility of several stellar populations in some of the 
sampled clusters from the presence of a c olor-extended 
HB, as suggested in some recent studies (D'Anton a fc | 
|Caloi 12008; Gratton et al.|201 0; Caloi & D'Antona 20fW 

A detailed explanation of the dereddening method is 
provided in Paper I, but we summarize here the differ- 
ent steps that we follow to apply the method, along with 
the results from every step, to obtain the final extinction 
maps and dereddened CMDs. The technique is applied 
independently to the two available colors {B — V and 
V — /) , and it is only at the end of the process that the 
extinction maps obtained from both analyses are aver- 
aged. 

The first step of our dereddening method is to calcu- 
late the probabilities for the observed stars to belong to 
the cluster. In order to do that, we first calculate the 
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surface density of the stars as functions of distance r to 
the cluster center , as explaine d in Paper I, and by a fit 
with an empirical King |l962) profile plus constant field 
model (see Figure [2] ana Ihble [6]) obtain, using Equation 
12 from Paper I, the conditionalprobabilities of the stars 
being members of the cluster as a function of position in 
the sky P{X = l|r) (see Figure [s]), where X indicates 
membership to the cluster (X = 1 indicates observa- 
tion of a member, and X = indicates observation of a 
non- member). Then we calculate, following Paper I, the 
probability densities of the observed stars as a function 
of r, color c, and magnitude m, and the probability den- 
sity of being field members as funct ion of c and m by 



comparison with a Besancon model ( [Robin et al. 12003 ) 
(see table [t]). Finally, applying Equation 20 from Paper 
I, we are able to find P{X = l|r, c, m), the conditional 
probabilities of the stars being members of the cluster as 
a function of position in the sky and in the CMD (see 
Figure [4] in this paper, and Figure 4 from Paper I as ex- 
amples of the method). We restrict our analysis to stars 
down to the completeness limit of our observations, and 
areas in which the probability P{X = l|r) > 0.1 (see Ta- 
ble |5|, as calculated in Equation 13 in Paper I. There are 
a few cases in which we do not calculate P{X = l|r, c, m) 
for one of the colors, since the observed CMD of the clus- 
ter is not calibrated in that color (NGC 6558 in 5 — V, 
and NGC 6235, NGC 6342 and NGC 6355 inV-I), and 
a comparison with the Besancon models can lead to mis- 
takes. In those cases, we equal the total probability to 
just P{X = l|r) (see Equation 13 in Paper I), and pro- 
ceed with the next steps of the analysis, which are not 
affected by any constant offset in the color calibration. 

Once the cluster membership probabilities have been 
calculated, we need to build the ridgeline for the CMDs 
of every cluster in our sample. We follow the three- step 
recipe described in Paper I to obtain the ridgeline^ for 
both available colors (see Figures [5]and|6|. 

A fter that, we move th e stars along the reddening vec- 



tor ([Schlegel et al.| T998| ) until they intersect the ridge- 
line, and smooth the resulting individual color excesses 
by building a bivariate non-parametric regression of the 
extinction as a function of the spatial coordinates right 
ascension and declination, which allows us to generate 
an extinction map|^ From this map, we take the relative 
extinction that corresponds to every star observed. We 
then plot the CMD of the stars in our observation after 
having been corrected for the differential reddening and 
use it as the input for the next iteration. 

Finally, after the process converges for both colors, we 
average the extinction maps obtained from both colors, 
generating a final extinction map and the dereddened 
CMDs (see Figures [t] to [3l]) . 

^ For some clusters the ridgelines were better calculated if we 
used only stars where the membership probability as a function of 
distance was P(X = l\r) > 0.5, instead of P(X = l|r) > 0.1 as 
in Paper I. For homegeneity on the process through the sample, 
we apply this more severe cut to calculate the ridgelines of all the 
observed clusters. 

^ For some clusters the reddening maps were better obtained 
when, in the iterative process described in section 2.4 of Paper I, 
we use stars with extinctions at most 2cr away from the extinction 
value calculated in the previous iteration, as opposed to the 3a 
threshold mentioned in Paper I. Again for homegeneity on the pro- 
cess through the sample, we apply this more severe cut to calculate 
the extinction maps of all the observed clusters. 



We now describe in more detail the individual cases 
of the different clusters in the sample. We display their 
B -V Ys. y and V - / vs. V CMDs, before and af- 
ter being differentially dereddened (see Figures [t] to [31]) , 
and describe any improvements in the photometry. The 
CMDs are plotted using the merged space and ground 
data, whenever space-based observations deeper than 
ground-based observations exist, or the ground photome- 
try in the rest of the cases. To quantify improvements in 
the photometry, we have measured the definition of the 
evolutionary sequences in the CMD by fitting Gaussian 
functions to different magnitude cuts in the different se- 
quences (1.5 mag above and 0.5 mag below the TO, this 
is, RGB and upper MS), similar to what we did in Fig- 
ure 9 of Paper I. To give an idea of the importance of the 
field star contamination, we compare the total number 
of observed stars in every FOV with the number of stars 
provided by the Besancon model of the Galaxy for the 
same pointing and FOV. 

We also present the extinction maps across the field 
of the clusters, along with their resolution and precision 
(see Figures [7|to|3l| and Table [8|. We compare these 
maps with tne SFD extinction maps (see Figure 32). 
The extinction maps provided by our dereddening tech- 
nique give the differential extinction in the field referred 
to the unknown absolute extinction where our ridgelines 
lie. The comparison with the SFD maps allows us to es- 
tablish these absolute extinction zero points. In order to 
make the comparison, we need to bring down our extinc- 
tion maps to the lower resolution of the SFD maps. Since 
our extinction maps are described by non-parametric 
continuous functions of the spatial coordinates right as- 
cension and declination, as mentioned in previous para- 
graphs, we evaluate these functions in ~ 1600 different 
points in every pixel of the SFD maps and average the 
obtained values using our precision maps to calculate a 
weighted average that provides us with one extinction 
value per SFD pixel. We observe that both our maps and 
the SFD maps generally agree in identifying the same re- 
gions of higher and lower extinction, although our maps 
usually show higher variation (reaching sometimes a fac- 
tor of 2) in the extinction values than the SFD maps. 
Notice also that the SFD maps give the extinction in- 
tegrated along the line of sight (i.e., extinction in the 
foreground, but also in the background of the observed 
clusters) , while the extinction maps provided by our tech- 
nique show the extinction just in the foreground of the 
observed clusters. 

4.1. NGC 6121 - M 4 (Figure^ 

The changes in the extinction across the field are sig- 
nificant {AE{B — V) = 0.15). The highest extinction is 
located in the northwestern region of the field, while the 
northeastern region of the field shows the lowest. There 
is a small blob of higher extinction {AE{B — V) ~ 0.05) 
~ 0.02 deg south of the center of the cluster, which 
the SFD map fails to show (see Figure 32). The ab- 
solute extinction zero point of our reddening map is 
E{B — V) = 0.50 from comparing our map with the SFD 
map. These significant values of extinction, both ab- 
solute and differential, are a consequence of the cluster 
being located behind the Ophiucus dust complex, about 
120 pc away from us. 

M 4 is likely the closest GC to the Sun. Because of this. 
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stars in the upper RGB are saturated in our images, even 
in the shortest exposures, and they do not appear in our 
CMDs of this cluster. The region studied for this clus- 
ter encompasses the whole FOV of our CCD, since the 
number of field stars present is small (< 20% of the stars 
in our CMD). Due to the proximity of M 4 to the Sun, 
most of the field stars are located behind the cluster. 
They are generally dimmer than the cluster stars, and 
the evolutionary sequences of the field stars are easily 
differentiable from the cluster sequences. The definition 
of the different branches in the CMD of the cluster are 
improved after the dereddening process, especially in the 
V^V — I CMD. The evolutionary sequences get ~ 35% 
tighter there. A noticeable feature in the CMDs of this 
cluster is the presence of a bimodal HB, with both red 
and blue components, though the blue HB (BHB) is not 
very extended. This bi modality of the HB, noted also 



in previous studies (e.g., Catelan et al. ( p!998) , and ref- 



erences therein), suggests the presence ot several stellar 
populations in this cluster. 

Some recent ground- based optical photometric studies 
of this cluster include Rosenberg et al.| (|2QQQa|) in VI 



Mochejska et al.| (|2QQ2|) m U BV , and [Ande rson et al 



in BV . Recent deep optical HST CMDs ot this 
cluste r have been presented in U^ ^eV^ 
eTH] ([20041), m by 



Ferdman et al. 



m U336 v^555-^8 i 4 by Bassa 
Richer et atj (| 2004 | , and 



and Hansen et al. 



( 2004 ) . This cluster is also a member of th e ACS survey 



of GGCs obtained in the V606^8i4 bands (Anderson et 
2008a) . R ecent infrared pho tometric studie s of this 



al 



cluste r are h y Cho fc Leej ( 2002]) in JK, a nd by 



Ferraro 



et al I ( |2000D in VJK. [Pulone et al.| ( |1999| also used the 
Ht>T for an optical ana infrared study of the stellar mass 
function of this cluster in V555/814 Jiio^ieo- The pres- 
ence of significant differential reddening in this cluster 
has been pointed out in many studies, an d a map of the 



differe ntial extinction was provided by , Mochejska et al. 
( 20021 ) ( see Figure p3v . Recently, the presence of distinct 



populations has been suggested from spectroscopic stud- 
ies of s tars in the RGB (Marino et al. 2008; Carretta et 
aIl [2Q10| ) and in the HB ( | Marino et al.|pDTTai r 

4.2. NGC 6I44 (Figure\^ 

NGC 6144 is also located behind the p Ophiuchi dust 
cloud and very close (- 40' northeast) to NGC 6121 (M 
4). Unsurprisingly, NGC 6144 also suffers from high and 
differential reddening. The variation of the extinction 
across the field is one of the highest shown in our sample 
{AE{B — V) = 0.5). Extinction is higher in the northern 
region of the observed field, ~ 0.1 deg from the center 
of the cluster, and peaks in the northeastern area. The 
absolute extinction zero point of our reddening map is 
E{B — V) = 0.75 from comparing our map with the SFD 
map. This is one of the few cases in which the differential 
extinction variation across the field is more important 
(by a factor of 2) in the SFD maps than in ours. 

The region studied for this cluster encompasses the 
whole FOV of our CCD, since there is only a moder- 
ate field star contamination (~ 50% of the stars in our 
CMD). Also the field stars are redder on average than 
the stars in the cluster whose different evolutionary se- 
quences are clearly di stinguishable in the CMD. Even if 
previo us works (e.g., "Martin & Whittet (1990); Green 



et al. ( 1992J ) give a higher value for the reddening vec 



tor towards p Ophiuchi dust cloud, with Ry ~ 4.3, still 
with the lower value assumed in our method, Ry ~ 3.1, 
the different evolutionary sequences get better defined 
(~ 25% tighter) in the CMD of the cluster after apply- 
ing the dereddening process. This cluster shows a short 
BHB, which suggests the presence of a single population, 
or the dominance of one population. 

The most recent ground-based optical photometric 
study of thi s cluster we found was done in BVI by |Neely| 
et al.) ( |2000 ') for the RGB and HB stars. From iJ gTdata, 
the first C MD of this cluster was published by ^Sa raje-| 
dini et al. ( 20Q7|) as a part of their 1^606^814 ACS sur- 
vey of GGCs. |Lan et al. ( 2010i) provide also CMDs in 
%36V606^<^656^675^814 using HST data. Infrared upper 
RGB p hotometry has been pres ented for this cluster by 
Davidge (2000) in JHK, and bylMinniti et al.|(|1995D in 
TTT 

4.3. NGC 6218 - M 12 (Figure^ 

The differential extinction across the field of this clus- 
ter is very mild {/S.E{B — V) ^ 0.04). Extinction is higher 
in the northern region (~ 0.1 deg of the cluster center). 
The absolute extinction zero point of our reddening map 
is E{B — V) = 0.18 from comparing with the SFD map. 

The region studied for this cluster encompasses the 
whole FOV of our CCD, since the number of field stars 
present is small (< 10% of the stars in our CMD). The 
mild differential extinction produces only a small im- 
provement in the CMD of the cluster after being dered- 
dened. The evolutionary sequences get ~ 10% tighter in 
y, y - /, a little less inV.B- V. Some of the CMD fea- 
tures, like the RGB (Thompson) bump, are more clearly 
visible now. A significant feature of this cluster is its 
BHB with a blue tail, which suggests the presence of 
several populations of stars. 

Ground-based VI photometric data on M12 down to a 
few magnitu des under the TO have most recently been 

fesented by von Braun et al. ( | 200 2 ) and R osenberg et| 
(2000b), B V data by [Brocato et a l. (1996),"WT^ata 
Hargis et al. (2004) and deeper V R data by de Marchi 
[eTal. ( 200 6|~g5T:B48 9 data on M12 were published 



Piotto et al. 



mm using the WFPC2. This cluster 
i s also a member of the ACS survey of GGCs in l^06^8i4 
( Anderson et al.|2008a ). We hav e found no infrared pho- 
tometrical study of this cluster, von Braun et al. (2002) 
mapped the extinction in this location, finding little dif- 
ferential reddening in the field (see Figure 33 ) . Recently, 
Carretta et al. (2007) have suggested the presence of two 
distinct populations with different helium co ntent from a 
spectroscopic study of RGB stars. Carretta et al. (2010) 
also find evidence for different populations from a spec- 
troscopic study of stars in the RGB. 



4.4. NGC 6235 (Figure\lL 

The changes in extinction across the field are moder- 
ate {AE{B — V) ~ 0.1), being the southeastern area 
of the field the region of highest extinction. The ab- 
solute extinction zero point of our reddening map is 
E{B — V) = 0.39 from comparing our map with the 
SFD map. These moderately high values of extinction, 
both absolute and differential, are a consequence of the 
cluster's projected position, close to the Ophiucus dust 
complex. 
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The region studied for this cluster includes only stars 
located less than 5.37' away from the cluster center, due 
to the high number of field stars relative to the GC stars 
(~ 85% of the stars when we look at the whole FOV of 
IMACS). As with NGC 6144, the field stars are redder 
on average than the cluster stars. This allow us to easily 
differentiate the cluster and field evolutionary sequences. 
These sequences in the cluster get ~ 15% tighter after be- 
ing differentially dereddened. This cluster shows a BHB, 
not very extended, which suggests the presence of a single 
population, or the dominance of one population. 

The most recent ground-based optical photometric 
study of this clu ster we found was don e in BV for RGB 
and HB stars by Rowland et al. (2003), who found a dif- 
ferential extinction relation along the X axis of their data 
and used it to correct their photometry. This cluster has 
also been observed with the WFPC2 aboard the iJSTand 
its ^439^^55 CMP has been published in the GGC survey 
bylPiotto et al." 



(2002|. Infrared upper RGB photometry 



has been presented tor this cluster by IDavicLge, (|^UUUj) m 
JHK and by Min niti et aT] ( |1995[ ) in JK. Deepe TJHK 
RGB near-i nfrared photometry was presented by |Chun 
eFaIl(l20T0l). 



4.5. NGC 6254 - M 10 (Figure 11) 



The extinction variation across the field of this cluster 
is mild {/S.E{B — V) ~ 0.1). Extinction is higher towards 
the west, peaking close ('^ 0.02 deg) to the cluster center. 
The absolute extinction zero point of our reddening map 
is E{B — V) = 0.28 from comparing our map with the 
SFD map. 

The region studied for this cluster encompasses the 
whole FOV of our CCD, since the number of field stars 
present is small (< 10% of the stars in our CMD). The 
CMD of this cluster, with an extended BHB, and its 
distance to the center of the Galaxy are very similar to 
those of M 12. The evolutionary sequences get ~ 25% 
tighter after being differentially dereddened, and features 
like the the RGB bump are more clearly visible now. THe 
BHB with a blue tail suggests the presence of several 
stellar populations. 

NGC 6254 has pr eviously been stud i ed usin g gro und- 
based faciliti es by Rosenberg et al. (|2000b|) and von 

by [Pollard et'at] (^200?) 



Braun et al.] (| 2002 |) m \/ 7 and 



BVl. [Fiotto &: Zoccalil ( |1999| ) used both ground- based 
facilities and Hbl ' for their study of this cluster in the VI 
bands. Beccari et al. ( ,2010) used also the HST for their 
mass segregation study of this cluster in the l^06^8i4 
bands. This cluster is also a member of the A C S survey 



of_GGCs_ji]_V606^8i4 (Anderson et al. 2008ak. Valenti 



a|. 

lira 



et al.| (l2004aj) studied this cluster in the infrared JK 
bands. I'he presence of differential reddening has been 
previously noted, and vo n Braun et al.| ([2002 ) mapped 
the extinction in the cluster field, finding a hig h er differ- 
ential extinction than in M 12 (see Figure [33| . |Carretta 
et al. (2010j) find evidence for different populations from 
a spectroscopic study of stars in the RGB. 



4.6. NGC 6266 - M 62 (Figure\R 

A band of material causing a significant increase in ex- 
tinction {/S.E{B — V) ~ 0.25) crosses the field in the east- 
west direction less than 0.02 deg south from the cluster 
center. The absolute extinction zero point of our redden- 
ing map is E{B — V) = 0.47 from comparing our map 



with the SFD map. In the SFD map the presence of this 
band of extinction is not so clear. 

The area studied for this cluster is restricted to the 
inner 6.84' of the cluster, due to the number of field stars 
(~ 50% of the stars when we look at the whole FOV 
of IMACS) and the relatively small tidal radius (r^ = 
11.28'). The definition of the different branches in the 
CMD of the cluster is highly improved, approximately by 
a factor of 2, after the dereddening process. Also now the 
presence of an extended blue horizontal branch (EBHB) 
is clearly visible and the location of the RGB bump is 
more easily identifiable. A striking feature in the CMD 
of this cluster is the presence of an extended and well 
populated HB, with both red and blue components. The 
BHB in particular is very extended. The shape of the HB 
and the fact that this is one of the most massive clusters 
in the Milky Way make this cluster a clear candidate to 
have several stellar populations. 

Previous work using ground-based data includes that 
of Rosenberg et al. (2000a) in VI, and Gerashchenk o & 

tlla (2004), Contreras et al. (2005) and Contreras e t al. 
\^'m BV. Be ccari et al.| ( |2006| ) studied this cluster 
using a combina tion oi BV 1 (ground) and /7255/7336V555 
(space) bands. Cocozza et al. (2008| studied a pulsar 
in M 62 and used their HST data to show the B — R 
vs. R CMD of the cluster. NGC 6266 is also included 
in the HST GGC s survey in publishe d by [Pi- 



rn me J^O- L o^cry s survey m ^3439 y/ 555 puunsne a uy ijri- 
otto et al.l (|2002D. | Valenti et al.| (|2007|) and Chun e^al. 



studied this cluster in the infrared JHK bands!" 
he presence of differential reddening has be en previ- 
ously noted, and Gerashchenko & Kadla ( 2004[ ) mapped 
the extinction towards the central region of the cluster 
(see Figure 33). 



4.7. NGC 6273 - M 19 (Figure\n 

The changes in the extinction across the field are sig- 
nificant {AE{B — V) ^ 0.3). Extinction is higher to- 
wards the east, and peaks in a small region located very 
close, less than 0.02 deg, of the center of the cluster. The 
absolute extinction zero point of our reddening map is 
E{B — V) = 0.32 from comparing our map with the SFD 
map. 

The area studied for this cluster is restricted to the 
inner 8.26' of the cluster, due to the number of field 
stars (~ 40% of the stars when we look at the whole 
FOV of IMACS). The definition of the different branches 
in the CMD of the cluster is highly improved after the 
dereddening process. The different sequences get tighter 
by approximately a factor of 2 after being differentially 
dereddened, and now the presence of an EBHB is clearly 
visible and the location of the RGB bump is more easily 
identifiable. The EBHB suggests the presence of several 
stellar populations in this cluster. 

The most recent optica l photometric study of this clus- 
ter we found was done by Piotto et al. (1999) in ^439^555 
from iJST images. As part of their study, tney obtained 
a map of the di fferential reddening of the central region 
of the cluster. Piotto et al." (2002") included the CMD 
from this work in their HST GGCs survey in ^43 9^555. 
In the infrared this cluster was studied in JHK by|Chun 
.et al. ( ,2010, ), by ,Davidge) ( ^2001^ ) and by Valenti et a 



4.8. NGC 6287 (Figure\14\ 
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The variation of the extinction across the field is one 
of the highest shown in our sample {AE{B — V) = 0.75). 
Extinction is higher in the northern region. A narrow 
band of material causing the highest extinction in the 
field crosses it in the east-west direction less than 0.02 deg 
from the cluster center. The absolute extinction zero 
point of our reddening map is E{B — V) = 0.67 from 
comparing our map with the SFD map. These significant 
values of extinction, both absolute and differential, are a 
consequence of the cluster's projected proximity to the 
Ophiucus dust complex. 

The area studied for this cluster is restricted to the 
inner 4.13' of the cluster due to the number of field stars 
75% of the stars when we look at the whole FOV of 
IMACS). The definition of the different branches in the 
CMD of the cluster is highly improved, approximately 
by a factor of 2, after the dereddening process. Now the 
RGB and SGB are easily identifiable, which was not the 
case before applying the dereddening technique. A short 
BHB is observed in the CMD, suggesting the presence of 
a single population, or the dominance of one population. 

Previous photometr ic stud i es in t he optical were done 
in BV by Stetson & West ( 1994|) from ground-based 
observations and in V 1 by .built on et al.| ( |1999| from 
space-based observations. NG(J 6287 is inc luded in the 
HST GGCs survey in ^439V555 publi s hed b y jPiotto et alT 



2002^ . In the infrared, Chun et al.|(|2010 |) and|Davidje 



^ ^ studied this cluster in J HK^ am Eee et sT. ( |20(]T] [ 

in Jiioi^ieo using the HST. 



4.9. NGC 6304 (Figure\l^ 

The variation of the extinction across the field is mod- 
erate {/^E{B -V)= 0.15). The area located - 0.05 deg 
east of the center of the cluster is the region of high- 
est extinction. The absolute extinction zero point of our 
reddening map is E{B — V) = 0.49 from comparing our 
map with the SFD map. 

The area studied for this cluster only includes stars 
with distances less than 3.65' away from the cluster cen- 
ter, due to the high number of field stars relative to the 
GC stars (> 90% of the stars when we look at the whole 
FOV of IMACS). The improvement of the CMD of this 
cluster is only marginal. The CMD of this cluster shows 
a very short, red HB (RHB), typical of high-metallicity 
GGCs. 

This cluster has been previously observed and studied 
using fflo und-based facilities by Qrtolani et al. (2000) in 
BV and Rosenberg et al. ( 2000a[ ) in VI. TTis one of 
the clusters in the GGCs survey by Piotto et al. (2002) 
i n ^489V555 and in the ACS GGCs survey in 1 /606-^814 



([And erson et al.||2008a[ ) done using HST data. "Valenti" 
et al. (2005) studied t his cluster in the infrared J HK 
EanHs. IValenti et al | ( [2QQ7|) i 
study in their infrared compilat 



include results from this 
compilation of bulge GGCs. 

4.10. NGC 6333 - M 9 (Figure\l^ 

The changes in the extinction across the field are signif- 
icant {/\E{B - V) - 0.25). The area located - 0.1 deg 
southwest of the cluster center is the region of highest 
extinction. The absolute extinction zero point of our 
reddening map is E{B — V)= 0.43 from comparing our 
map with the SFD map. 

The area studied for this cluster is restricted to the 
inner 5.59' of the cluster due to the number of field stars 



(~ 45% of the stars when we look at the whole FOV of 
IMACS) and the relatively small tidal radius (r^ = 8.16'). 
The evolutionary sequences of M 9 get ~ 35% tighter af- 
ter being differentially dereddened. This cluster posseses 
a BHB with a blue tail, which suggests the presence of 
several stellar populations in this cluster. 

The most recent optical photo metric study o f this clus- 
ter we found is that of Janes & Heasley ( 1991[ ) in BV for 
stars brighter than the I'U point. Although this cluster 
has been previously observed with the HST (see Table 

t, the first CMD of M 9 using space data is published in 
is paper in BVI . In the infrared jDavidge| (|200l|) and 
Chun et al. (2010) studied this cluster in J J 



4.11. NGC 6342 (Figure 17) 



The changes in the extinction across the field are sig- 
nificant {AE{B — V) ~ 0.4). The northwestern quadrant 
of the field, reaching the projected position of the center 
of the cluster, is the region of highest extinction. The 
absolute extinction zero point of our reddening map is 
E{B — V) = 0.59 from comparing our map with the SFD 
map. 

The area studied for this cluster reaches only stars with 
distances less than 4.21' away from the cluster center, 
due to the significant number of field stars relative to 
GC stars (~ 70% of the stars when we look at the whole 
FOV of IMACS). The definition of the different branches 
in the CMD of the cluster is highly improved after the 
dereddening process, getting tighter by a factor of 2. Now 
the RGB and SGB are easily identifiable, along with a 
RHB, which was not the case before applying the dered- 
dening technique. The RHB is typical of high-metallicity 
GGCs. 

The only recent optic al photometric studies fou nd for 



this cluster are those of Heitsch & Richtler 



1999) in VI 



lotto et al. 



using ground-based facilities and that of 
(2002) in ^439^555 using the HST. In the infrared Val enti 
et al. (2004b) studied this cluster in VJHK. ValentTet 
al. (200^ include results from this study in their infrared 
compilation of bulge GGCs. A map of the differential 
reddening was produced by Heitsch & Richtler (1999) 
(see Figure 33). 



4.12. NGC 6352 (Figure\l^ 

The differential extinction variation across the field of 
this cluster is only moderate {AE{B — V) ^ 0.15). Ex- 
tinction is higher in the northern region of the analyzed 
field. The absolute extinction zero point of our redden- 
ing map is E{B — V) = 0.34 from comparing our map 
with the SFD map. 

The area studied for this cluster is restricted to the 
inner 5.61' of the cluster, due to the high number of 
field stars relative to the GC stars (> 90% of the stars 
when we look at all the available FOV). The evolution- 
ary sequences get ~ 30% tighter after being differentially 
dereddened. The cluster shows a RHB, typical of high- 
metallicity GGCs. 
This clus ter was studied by Saraiedini fc Norris ( 1994|) 



in BV and Rosenberg et al. 



( |2000ap and Fancino et aL 
d, and from space byll^ull- 



ton et al. 


(1995|) in 


by 'l^ aria ( 


l^eltzing 



the ACS 
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survey of GGCs in Feoe^SM C Ande rson et aL|2QQ8aD . We 
hav e found no infr ared photometrical study of this clus- 
ter. i Pancino et al.| (|2010j) show the presence of two pop- 
ulations of stars with anticorrelated CN and CH band 
strengths. 

4.13. NGC 6355 (Figure\l^ 

The changes in the extinction across the field are sig- 
nificant {AE{B — V) ^ 0.3). The southwestern area 
of the cluster is the region of highest extinction. The 
absolute extinction zero point of our reddening map is 
E{B — V) = 1.20 from comparing our map with the SFD 
map. 

The area studied for this cluster is restricted to the 
inner 2.36' of the cluster, due to the high number of 
field stars relative to the GC stars (> 90% of the stars 
when we look at all the available FOV). The evolution- 
ary sequences get ~ 10% tighter after being differentially 
dereddened. This cluster posseses a BHB, although the 
small number of stars in the HB, due to the small stud- 
ied field, makes it difficult to infer any information about 
different populations in this cluster. 

The only previo us recent optica l photo metric studies 
found are those by Ortolani et al. ( 2Q03|) in BVI us'm^ 
ground-based facilities, and that by Fiotto et al.' ('2002) 
in B/iRci y^^^ using the HST. In the infrared Chun et al. 



i n i343 C) \/555 usm^ une no i. iii uiie iiiirarea ^^iiuii eu ai. 
( |20 10| , |Davidge| (|200 1 [ ) , and |Valenti et al.| ( |2007| ) studied 
this cluster in JH K . 



4.14. NGC 6397 (Figure\2§ 

The differential extinction in this cluster is only mild 
{AE{B — V) ~ 0.05). Extinction is not especially concen- 
trated in any area of the field. The absolute extinction 
zero point of our reddening map is E{B — V) = 0.18 from 
comparing our map with the SFD map. 

This is the cluster with the smallest apparent distance 
modulus, though it is not the one closest to the Sun. As 
for M 4, stars in the upper RGB are saturated in our 
images, even in the shortest exposures, and they do not 
appear in our CMDs of this cluster. The contamination 
from field stars is moderate (~ 30% of the stars when we 
look at all the available FOV) , mainly affecting MS stars 
in the cluster. This moderate contamination allows us to 
study almost all the available FOV, up to 8.18' from the 
cluster center. The sequences of NGC 6397 get ~ 10% 
tighter after being differentially dereddened. This cluster 
shows a short BHB, which suggests the presence of a 
single population, or the dominance of one population. 

Some of the most recent studies of thi s cluster done 
with ground-base d fac ilities are those by [Rosenberg et 



al. (200 0a|) in 1//, [Anderson et al. (2006| in UBVI, and 
n^i„^^^. ^1 I /lony^p BT/T^rhis cluster has also been 
1TST. A series of papers has been pub- 



Kaluzny et al. 
observed with t' 
lished dedicated to the study of the characteristics of this 
cluste r using very dee p phot ometry obtained with the 
ACS ([Ande r son et al . 2008b; Hurley et al.[[2008| [Richer 
^!008l Ipavis et al. 2008). NGC 6397 is also^n- 
lie GGC survey from Fiotto et al. (2002)in 
and of the ACS GGC survey in 1/69 6-^314 ( |An- 



et al. 

eluded in t 



derson et al. 2008a). 



(20001 also used 



de Marchi et al. 
the Ht>T for an optical and infrared study ot' the stellar 
mass function of this cluster in V606^8i4-^iio^i60- Recent 
spectroscopic studies of this cluster have found different 



chemistry among its stars: Carretta et al. (2010) show 
some 0-Na anti-correlation is present, although the num- 
ber statistics is small (onl y 16 stars studied, t he smallest 
in their sample of GCs); [Lind et al. ( [201 1 ) claim the 
presence of two different stellar populations that have 
different chemical compositions of N, O, Na, Mg, and 
probably Al, and note that the cluster is clearly domi- 
nated (75%) by the second generation. 



4.15. NGC 6522 (Figure^ 

The differential extinction variation in this cluster is 
moderate {/S.E{B — V) = 0.1), although we should realize 
that the studied region is also small (see next paragraph). 
Extinction is higher in the southwestern region of the 
analyzed field. The absolute extinction zero point of our 
reddening map is E{B — V) = 0.58 from comparing our 
map with SFD map. 

This cluster is located in Baade's Window. The high 
number of non-cluster member stars (> 90% of the stars 
when we look at all the available FOV) restricts the study 
in this region to stars closer than 2.33' from the cluster 
center. The improvement of the CMD of this cluster after 
being dereddened is only marginal. It shows a BHB with 
a blue tail, which suggests the presence of several stellar 
populations in this cluster. 

The most recent optical photometric s tudies found, 
done with ground-based da ta, are those by[Barbuy et al.[ 
^94) in BV, and Terndru pfc Walker[ ( |l994p m BVI 
[T^rndrup et al. (1998) show a differential extinction rela- 
tion with respect to the center of the cluster and use it to 
correct their previous phot ometry . Usin g the HST^ the 
cluster was studied byjSh^xaetah ( 1998 ) in the U^^qB/^2,9 
bands usin g pre-COS'lAK data .Tt is also one of the clus- 
ters in the [Fiotto et al.' (2002) GCs survey in ^439^ 555 - 
In the infrared, [Pavidg e (2001) and [Valenti"et al.,(201Q| 
studied this cluster m JHK~ 



4.16. NGC 6541 (Figure^ 

The differential extinction in this cluster is moderate 
{AE{B — V) ^ 0.1). Extinction is not especially con- 
centrated in any area of the field, although there seems 
to be an increase towards the southeast area of the ob- 
served field. The absolute extinction zero point of our 
reddening map is E{B — V) = 0.16 from comparing our 
map with the SFD map. 

The region studied for this cluster encompasses most of 
the FOV of our CCD, up to 8.63' from the cluster center, 
since there is only a moderate field star contamination 
(~ 40% of the stars in our CMD). Also the field stars are 
redder on average than the stars in the cluster, which 
lets us distinguish NGC 6541 evolutionary sequences un- 
ambiguosly. The sequences in the cluster get ~ 10% 
tighter after being differentially dereddened. This clus- 
ter possesses a BHB with a blue tail, which suggests the 
presence of several stellar populations. 

Ground- based s tudies of this cluster w ere presented 



( [1997| ) in BVRI and by [Rosenberg et 
It has al so been studied using Htl'I 



by Alcaino et al. 
al.| (|2000aD ^mT 

WFFC2 data by [Lee fc Ca xney (2006) in VI, and it is 
one of the clus ters in th e ACS GGC survey observ ed in 

In the infrared 
in JHK. 



V6n6-^8]4 (Anderson et al.^^2008a) . 
et al. ( [2006p studied this cluster' ir 



Kim 



4.17. NGC 6553 (Figure 23) 
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There are significant changes in the extinction across 
the field {AE{B - V) - 0.25). In our maps, the high- 
est extinction is located in the northern region of the 
field, while the southwestern region of the field shows 
the lowest, which is not the trend shown in the SFD 
maps, but it does a^ee well with the trend of the map 
provided by [Heitsch fc Richtler] ( |1999| ) (see Figure [33| . 



The absolute extmction zero pomt of our reddening map 
is E{B — V) = 1.35 from comparing our map with the 
SFD map. 

The area studied for this cluster is restricted to the 
inner 3.91' of the cluster, due to the high number of field 
stars relative to GC stars (> 90% of the stars when we 
look at all the available FOV). The upper MS and RGB 
get ~ 30% tighter after being differentially dereddened, 
and the RGB bump now is more clearly identifiable. The 
well-populated bump, the very extended upper RGB and 
a RHB, suggest that this is a metal-rich GGCs. 

Some of the most recent optical photometric studies 
of this cluster are that by [Sagar et al.[ (|1999|) in VI 
grou nd-based data and t hose by peauliei 



usm^ 



data. 



(2001) and ZoccaU 



ZoccalTet al 



ii et al I ( [200 
. ( 2dnTf calci 



leu et al. 

2m) in F555/814 using HST 



culated the extinction map, 
not sh own in their paper , and applied it to get a clean 



Guarnieri et al. 



CMD. Heitsch fc Richtler ( 1999j) show an ex tinction map 

S see F ig ure [33|) calculatea with data from Sagar et al. 
19991. 



( Il998| ) in VIJK, and [l^'erraro 
letaJj (200(^ m VJ K studied this clust er in the optical 
and in The infrared. Valenti et al.| ( [2QQ7 ) includes results 
from this last study in their infrared compilation of bulge 
GGCs. 

4.18. NGC 6558 (Figure^ 

The extinction variations in this cluster are significant 
{AE{B — V) ^ 0.15), especially when we consider the 
small size of the studied field (see next paragraph). The 
absolute extinction zero point of our reddening map is 
E{B — V) = 0.51 from comparing our map with SFD 
map. 

The high number of non-cluster member stars (> 90% 
of the stars when we look at the whole FOV of IMACS) 
restrict the study of this region to stars closer than 1.81' 
from the cluster center. The improvement of the CMD 
is only marginal. This cluster shows a BHB, although 
the low definition of the HB due to the small number of 
stars in this region caused by the small field FOV studied 
makes it difficult to infer any information about different 
populations in this cluster. 
We found only two previous optical photometri c stud- 



ies usi ng ground-based data , those by Rich et al. ([1998) 
and by Barbuy et al] (|2007[), both in VI. Although this 



cluster has been previously observed with the HST^ the 
first CMD of the deep MS stars using space-ba sed data 
is publish ed in this paper in VI. In the infrared [Chun etj 
ad (I20T0I studied this cluster in JHK. 



4.19. NGC 6624 (Figure\25 

The differential extinction across the field is moderate 
{AE{B — V) ~ 0.15). Extinction is higher in the south- 
eastern region of the analyzed field. The absolute extinc- 
tion zero point of our reddening map is E{B — V) = 0.25 
from comparing our map with the SFD map. 

The area studied for this cluster is restricted to the 
inner 3.92' of the cluster, due to the high number of 



field stars relative to the GC stars (> 90% of the stars 
when we look at all the available FOV). There are also 
some dim blue field stars in our CMD that belong to 
the Sagittarius dwarf galaxy. Stars from this galaxy's 
MS are easily observed in the CMD at magnitudes dim- 
mer than V ~ 21 and colors bluer than the cluster MS. 
The improvement of the CMD of the cluster after being 
dereddened is only marginal. This cluster possesses a 
RHB, characteristic of a metal-rich GGC. 

This cluster has been s tudied using ground-based data 
by [Sarajedini fc Norris (1994), who presente d a photo- 



met ric stu dy o t the KGij and HB stars in BV. [Richtler et 



al. ( [1994 ) and Rosenberg et al. (|2000a') obtained deeper 
pnotometry, reaching a tew magnitudes below the TO 



in BV and VI ^ respectivel y. Using HST, studie s of this 
cluster were published by Sosin fc King [ ([1995 [) in BV 
using pre-COSTAR data, and by Guhathakurta et al 
(1996) in /73365439^555 and Heasley et al. ( 2000 ) in 



This cluster is also a member of the Fiotto et al. (2002) 
^489^55 5 GGCs survey, and of the V6n6-^8i4 A(J! S GGCs 



survey ( Anderson et al. 2008a). In the infrared Valenti 
'et al.l ( l^04b|) studied this cluster in VJHK. 'Valenti et 
^al., ( ^2007|) include results from this study in their infrared 
conipilation of bulge GGCs. 



4.20. NGC 6626 - M 28 (Figure 26) 



The extinction variations across the field of this cluster 
are significant {AE{B — V) ^ 0.2). Extinction is lower 
at the projected center and southwestern region of the 
observed field, which is not the trend observed in the 
SFD maps. The absolute extinction zero point of our 
reddening map is E{B — V) = 0.46 from comparing our 
map with the SFD map. 

The area studied for this cluster is restricted to the 
inner 4.95' of the cluster, due to the high number of field 
stars relative to the GC stars (> 90% of the stars when 
we look at all the available FOV). Due to its proximity, 
the cluster's brightest stars are brighter on average than 
the field's, and we can identify without much difficulty 
stars from the RGB. The evolutionary sequences of M28 
get ~ 35% tighter after being differentially dereddened. 
This cluster shows a BHB with a blue tail, which suggests 
the presence of several stellar populations. 

This cluster has be en studied in th e optica l using 
ground-based data by ^ose nber^ et al.J (2000a) in VI, 
and u s ing i/ grdata byp -esta et al. (200O and by Golden 
et al. (2001) in VI. [Davidge et al., ( |199d|) used a combi- 
nation of optical and infrare d bands {Bv JK) to study 
this cluster. In the infrared Chun et al. (2010) studied 
this cluster in JHK. 



4.21. NGC 6637 - M 69 (Figure^ 

The differential extinction across most of the observed 
field is moderate {AE{B — V) ^ 0.1). Extinction is not 
especially concentrated in any area of the field. The 
absolute extinction zero point of our reddening map is 
E{B — V) = 0.17 from comparing our map with the SFD 
map. 

The area studied for this cluster is restricted to the 
inner 4.39' of the cluster, due to the high number of 
field stars relative to the GC stars (> 90% of the stars 
when we look at all the available FOV). There are also 
some dim blue field stars in our CMD that belong to the 
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Sagittarius dwarf galaxy. Stars from this galaxy's MS 
are easily observed in the CMD at magnitudes dimmer 
than V ^ 21 and colors bluer than the cluster MS. The 
improvement of the CMD of this cluster is only marginal. 
This cluster posseses a RHB, typical of high-metallicity 
GGCs. 

This cluster has been the subject of several studies, 
usually along with the also metal-rich NG C 6624. Us- 
ing fflound-based data Ferraro et al. (1994), in BVJK, 



and Sarajedi ni fc Norrisl (|l994l ), in i^V, present a pho- 
tometric stud y ot the KG B an a HB stars of this cluster, 
and Richtler et al.| ( |1994l ) and [Rosenberg et al. (2000aj) 
got deeper photometry, reaching a few magnitudes below 
the TO, in BV and VI respectively. Using HST data, 
Heasley et al. ( 2000 ) studie d this cluster in VI- This clus- 
ter is also a rnember of thejPiotto et al. (2002]) ^439^555 
GGCs survey, a nd of the V606^8i4 ACS^GGCs survey 
( Anderson et al.|[2008aK In the infrared this clu ster was 
studied by" Ferraro et al.j ([2000| in VJK and by Valenti 



et al.| ( |2005|) m JHK . 'Valen ti et al.| ( |2007 !) include re- 
sults from these studies in their infrared coe 



bulge GGCs. 



compilation of 



4.22. NGC 6642 (Figure^ 



The differential extinction in this cluster is high 
{AE{B — V) ^ 0.2). The absolute extinction zero point 
of our reddening map is E{B — V) = 0.36 from comparing 
our map with the SFD map. 

During the observation of this cluster with the Magel- 
lan telescope, we had the least photometric conditions of 
our run (see Table 11). This resulted in the shallowest 
CMD of our study, reaching only stars down to the very 
upper parts of the MS. The high number of non-cluster 
member stars (> 90% of the stars when we look to all the 
available FOV) restricts the study in this area to stars 
closer than 2.22' from the cluster center. The improve- 
ment of the CMD of this cluster after being dereddened 
is only marginal. A prominent feature in the CMDs of 
this cluster is the presence of a bimodal HB, with both 
red and blue components, though the blue BHB is not 
very extended. This bimodality of the HB can suggest 
the presence of several stellar populations. 

The most recent optical photometric stud y of NGC 
6642 using ground-based data was done by Bar buy et al. 
([2006) in BVL Using i7gr, Balbinot et al (20Q9| stud- 
ied this cluster in 1 ^606-^81 4- iNGC 6642 is also a member 



of the[Piotto et al.'(^002) ^439^555 GGCs survey . In the 
infrareTpJ avid^e (200^7 Kim et al. (2006), and Valenti 
et al. ( 2(7(J7p studied tms cluster in JH K . 



4.23. NGC 6656 - M 22 (Figure^ 

The extinction variation across the field of this cluster 
is moderate {AE{B — V) ~ 0.1). Extinction is lower in 
the southern area of the field. The absolute extinction 
zero point of our reddening map is E{B — V) = 0.33 from 
comparing our map with the SFD map. 

The relative proximity of M22 explains why the upper 
RGB stars are saturated in our observations and are ex- 
cluded from the CMDs of this cluster. The region studied 
for this cluster encompasses the whole FOV of our CCD, 
since there is only a moderate field star contamination 
40% of the stars in our CMD). Due to the proxim- 
ity of M 22, most of the field stars are located behind 



the cluster. They are generally dimmer than the cluster 
stars, and the evolutionary sequences of the field stars 
are easily differentiable from the brighter sequences of 
the cluster. After being dereddened, the sequences of 
the cluster get ~ 25% tighter in V, V — /, a little less 
in V^B — V . An EBHB suggests the presence of several 
stellar populations in this cluster. 

There are several optical photometric studies for this 
cluster. Using ground-based data, it was studied by 
120041) in By/ bylWebb et al. 



Monaco et al. 



UV.by Kaluzny ^ 'Fhompson 



Tl 



berg et aE 12000a) in , and b y Anthony-T warog e t al. 
(|l995[), by |Richter et al.| ( |1999[ ) and bylLee et al. (20^^ 



by iWe bb et al.i (|2 004|) in 
( MOT| inijl/,i )y|Ro^^^^^ 



in Str omgrem uvbyCa, vby a nd VbyCa bands, respec- 
tively. Piotto & Zoccali (1999) used both ground facilies 
and the Ht) T for their study of this cluster in V/, and 
Anderson et al. ( 2003| ) show CMDs for the cluster in 
^435^606^675^814 bands. This cluster is also a member 
of GGCs in ^606^814 



Anderson et al. 



of the AC S survey o f GGCs in ^606^814 (1^ 

[2Q08a| ). iPiottol d2QQ9| ) shows a double SG B m his stuHy 
of M 2 2 using 1^606-^81 4- I n the infr ared [Minniti et al. 



( [T995I ), [Davidge fc Harris | (1996) and C ho fc Lee| ( |2002p 



studied this cluster in JK . The presence of a metallicity 
spread in this cluster has been the subject of a long de- 
bate, co mplicated by the pres ence of inster stellar redden- 
ing (see Pancino et al. (2010) for a brief historic review). 
Some recent studies seem to confirm the presence of at 
least two stellar populations with different chemical com- 
positions, both photometrically (Piotto 200 9| and spec- 
troscopically ( Marino et al. 2009; Da C osta et al.||2009" 



i al.||i : 

studie s by [Kayser et al., (^2008) and by ^ 

CH and CN anti-correlations or 



Pancino et al. 



(2010) do not tind clear 



Carretta et al.||2010; Marino et al. 2011b)7^ItEougFtl£ 



bimodalities for the stars in their sample. 

4.24. NGC 6681 - M 70 (Figure^ 

The differential extinction across most of the observed 
field is mild {AE{B - V) - 0.05). Extinction is not 
especially concentrated in any area of the field. The 
absolute extinction zero point of our reddening map is 
E(yB — V) = 0.09 from comparing our map with the SFD 
map. 

The area studied for this cluster is restricted to the 
inner 5.21' of the cluster, due to the high number of 
Galactic field stars relative to the GC stars (~ 70% of the 
stars when we look to all the available FOV). The field 
stars are redder on average than the cluster stars, which 
allows us to easily distinguish the different evolutionary 
sequences. There are also some dim blue field stars in 
our CMD that belong to the Sagittarius dwarf galaxy. 
Stars from this galaxy's MS are easily observed in the 
CMD at magnitudes dimmer than V ~ 21 and colors 
bluer than the cluster MS. The improvement of the CMD 
of this cluster is only marginal. M 70 posseses a BHB 
with a blue tail, suggesting the presence of several stellar 
populations in this cluster. 

The most recent optical photometric studies we found, 
using ground-based d ata, a re those by Brocato et_aL] 
(1996) in BV, a nd by [Rosenberg et ah (2000a) in 



Using iJgr data, [Piotto et aL| ( |2002p included it in their 
GGCs survey in i:>435 1/555 • This cluster is also a mem- 
ber of the V606^8i4 ACS survey of GGCs ( Anderson et 



ar|[2008a[ ). In the infrared [Kim et al.| ( [2006[ ) studied this 
cluster m JHK. 
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4.25. NGC 6809 - M55 (Figure^ 

The extinction variation across the field of this cluster 
is only mild {AE{B — V) ^ 0.05). Extinction is higher 
in the northwestern area of the observed field. The im- 
provement of the CMD of this cluster is only marginal. 
The absolute extinction zero point of our reddening map 
is E{B — V) = 0.14 from comparing our map with the 
SFD map. 

The region studied for this cluster encompasses the 
whole FOV of our CCD, since the number of field stars 
present is small (< 10% of the stars in our CMD), though 
there is clear contamination from stars in the background 
Sagittarius dwarf galaxy. Stars from its MS are easily ob- 
served in the CMD at magnitudes dimmer than V ~ 21 
and colors bluer than the cluster MS. The tip of the RGB 
is saturated in our CMDs of this cluster, due to the small 
apparent distance modulus of M 55. The improvement 
of the CMD of this cluster is only marginal. The cluster 
shows a BHB, not very extended, suggesting the presence 
of a single population, or the dominance of one popula- 
tion. 

The latest optical photometric studies of this clus- 
ter done using ground-based data are by JKaluzny et al.| 
([2QQ5| in UB VI, by |Pych et ar|(|2001l) andby Zloczewskl 
^^^^ ([2011") in BV, and by |KQsenberg~et al.| (|2000a^ 
. Using; a combination of around- based and Ht)T 



et al. 

data, this clu ster was studied by Lanzoni et al. (12007) in 



BVI, and by Bassa et ah 



i20mm~BVRlM r.. the ACS 



VI photometry used by bassa et al. (2008) comes from 



Anderson et al. (2008a), since this cluster is one of the 
members of the ACS GGC survey. In the infrared this 
cluster was studied by Ferraro et al. (2000) in VJK. The 
presence of m ore than one s t ellar population has been 
suggested by |Carretta et al. (2010) based on the Na-0 
anti-correlation fou nd in th eir study, but the stu dies by 



Kayser et al | ([20081 ) and by [Pancino et al.| ( [2Q1Q| ) do not 
find a clear CH and CN anti-correlations or bimodalities 
on the stars in their sample. 



5. SUMMARY 



This article is the second of a series devoted to the 
study of the inner GGCs. High and patchy differential 
extintion in the line of sight makes difficult to extract 
accurate and meaningful information about the physi- 
cal parameters of the GCs located in the Galactic bulge 
and inner disk and halo. In this article we present a 
highly homogeneous and complete optical sample of the 
brightest GCs located towards the center of the Milky 
Way. We produce some of the deepest CMDs of the 25 
inner GGCs that comprise our sample, using Magellan 
6.5m telescope and HST observations. We apply a new 
dereddening technique to map the differential extinction 
in front of these GCs and eliminate their effects in the 
CMD0 The GGCs in our sample show tighter evolu- 
tionary sequences, up to a factor of 2 in some cases, in 
the CMDs built after applying the dereddening method. 
These cleaned CMDs will be used in subsequent papers 
to analyze the astrophysical parameters of the GCs and 
their stellar populations. The high-definition extinction 
maps will be also used in future papers to study the prop- 
erties of the obscuring material along the GGCs' lines of 
sight. 
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Institute of Technology, funded by the National Aero- 
nautics and Space Administration and the National Sci- 
ence Foundation. This research used the facilities of the 
Canadian Astronomy Data Centre operated by the Na- 
tional Research Council of Canada with the support of 
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APPENDIX 
DOPHOT 



In this work, we have obtained the photometry using a new version of DoPHOT. DoPHOT is a software created to 
extract PSF photometry from astronomic al images. A com plete and deep explanation of its main characteristics, and 



a description of how it works, is given by Schechter et al. 



1993) 



The main new characteristics of the distribution of DoPHOT we developed for this work are: 

• All the program subroutines and functions have been rewritten to make them compatible with the standards of 
Fortran90. 



• We have tried to make the subroutines and functions more easily readable. To achieve this, we have eliminated 
all the common statements present in previous versions and named all the variables used by the subroutine and 
function in the call statement. Also with the same spirit of making the program more understandable we have 
created, combined or eliminated some of the subroutines, and included, at the beginning of each file, a brief 
explanation of what every function and subroutine does. 

• The program works now with real values for the pixel values, so there is no need to worry about conversions to 
integers as in previous distributions. 



^ We are working on making this photometric database easily 
accesible to the community through a webpage 
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• The stars are masked, obliterated, added, and subtracted across the different subroutines of the program in 
circular sections, and not in rectangular sections as before. This way fewer usable pixels are obliterated, and 
fewer pixels outside an interesting radius are dealt with. 

In addition to this, this distribution shares some of the characteristics of the program that were present in previous 
distributions, but that were not mentioned in the original paper: 

• DoPHOT allows for a variable PSF model in the image. 

• DoPHOT can perform an aperture correction to the magnitudes obtained, according to their positions and 
magnitudes. The program checks for variations in the aperture correction in magnitude and in position in the 
image, and corrects the fitted magnitude accordingly. 

This new version of DoPHOT is publicly ava ilable through: 



htt p : / / www . astro . puc . cl/ ^ j alonso / dophot .tar 



INSTRUMENTAL CALIBRATION 



The camera in IMACS consists of eight chips, each with its own slightly different mean quantum efficiency and 
color sensitivity. To calibrate our data, we must first transform the photometry from each of the individual chips t o a 
common zero point, then transform this system to a standard system-in our case, the Johnson-Cousins system (Be ssell| 
p^979) . In this appendix, we shall focus on the first step of transforming the individual chips to a common instrumental 
system. 



For this purpose, we arbitrarily chose one of the chips, chip 2 (see Figure 34), as the one defining the instrumental 
system. The goal was to tie all the other chips to this chip's photometric zero point. We attempted to do this in two 
ways. First, we observed one Landolt (1992) field (PG 1047+003) in every one of the chips and in all three filters (BVI) 
in rapid succession on every night of our run. There were very few stars in this field, and the scatter we found in the 
mean photometric offsets - about 0.05 mag between chips - was too large given our need to bring all our photometry 
to a common system to 1-2% precision. We surmise that the relatively poor results from this approach reflects the 
quality of the nights during the 2005 run, and also possible flat-field variations at the 2-3% level for this run. 

We tried a second approach during a sequence of observations during a 2006 Magellan run (see Table [2| in which we 
tried to auto-calibrate the IMACS chip array using stars from the cluster NGC 6397. The basic idea was to obtain a 
minimal number of exposures in each filter so that at least part of all chips overlapped with chip 2, the chip defining 
our instrumental photometric system. We found that a sequence of three exposures could achieve the desired overlap. 



Figure 34 illustrates the approach. In the left panel, the darker outline of the IMACS chip shows the location of a 
first exposure, in this case at the same location as the actual field we observed in 2005 for the cluster NGC 6397. The 
second exposure, shown as a lighter (red) outline is offset in both RA and Dec so that four of its chips (1, 2, 5 and 6) 
overlap with the position of chip 2 in the first exposure. The right panel again shows the initial pointing in the darker 
outline, with a third pointing as a lighter (green) outline. Now, chips 3, 4, 7 and 8 of the last exposure overlap with 
chip 2 from the first exposure. This sequence, repeated for each filter, generated a large list of stars that have been 
observed on chip2 and on every other chip in the array. Because only three exposures were required and because we 
took these in rapid succession, we can safely ignore airmass variations as we calculate the offsets from a given chip to 
the system defined by chip 2. In addition, we have exposures from the 2005 and 2006 runs of the same fields (the dark 
outlines in Figure [34|) that allow us to determine any possible variations in chip sensitivities between runs. 
The model we used to transform each chip to the chip 2 system has the form 

1712 = -\- Zi -\- JiCi^ (Bl) 

where rrii is the instrumental magnitude of chip i (excluding i = 2), is the zero-point offset from chip i to chip 
2, and ji is the color coefficient for instrumental color q. In practice, we carried out this analysis separately for 
each filter, B, V, and I, with the colors defined appropriately for each filter {{B — V) for the B-band transformation 
and {V — I) for the V and I-band transformations). We found that the color terms are all negligibly small in the 
chip-to-chip transformations for each filter, so the final model adopts 7i = for all chips and all colors. Thus we only 
had to detemine the mean values of the Zi coefficients as the weighted mean of the differences of the instrumental 
magnitudes, — m2, for all stars in common between chip i and chip 2. Weights for each measurement were assigned 
according to the photometric errors and the dispersion in the offset ditribution. A first application of this approach. 



to the calibration sequence of NGC 6397 obtained in 2006 (Figure 34), allowed us to place all the data from that run 
on the 2006 chip 2 instrumental system. 

Since the relative sensitivities of the chips might have changed between the 2005 and 2006 runs, we then took an 
additional step to define the chip-to-chip offsets for the 2005 run from which all our ground-based cluster photometry is 
defined. This process involves the observations of NGC 6397 from our 2006 run (described above) and the observations 
of the same field for the same cluster obtained in 2005. 

For this purpose we adopted the following transformation: 

^i,06 = ^i,05 + + ^iCi + koB^ob - koeXoG' (B2) 

The zero-point offset from the 2005 observations required to place them on the 2006 system for chip i is given by 
Yi (this now includes the case z = 2), for a color term J^, and an instrumental color on the 2005 system of q. The 
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final two terms account for the different airmasses and possibly different extinction coefficients for each run and each 
observation. We do not know this product (extinction coefficient k times airmass X) for the 2006 run, because we did 
not observe standard stars on that run. But we know that the sum of these products is constant when comparing the 
2005 and 2006 NGC 6397 observations. Thus, we can rewrite the equation above as 

m,06 = m,05 + + ^iCi, (B3) 
in close analogy to Equation Bl. Note that the new zero point offset now implicitly includes the airmass terms: 

= Y,^ ^05^05 - ^06^06. (B4) 

We can simplify this due to the fact that we found that the color terms in Equation B3 are all negligibly small. Thus we 
only had to detemine the mean values of the Y- coefficients as the weighted mean of the differences of the instrumental 
magnitudes, m^^oe — ^i,05 5 for all stars from chip i. As before, weights for each measurement were assigned according 
to the photometric errors and the dispersion in the offsets distribution. 

At this stage, we can transform the 2005 NGC 6397 data to a common instrumental system defined by chip 2 from 
the 2006 observations of this same cluster. But what we really require are the offsets in the 2005 system that bring all 
data from the IMACS chips obtained in 2005 to a common system. That is, the final chip-to-chip offsets for the 2005 
data are simply 

ay: = Y: - Y^. (B5) 

The values of Z and AY^ required to tie the various chips systems to a common system defined by chip 2 are listed 
m Table [9| The errors in the values of Z and AY"' are a ~ 0.001, the error of the mean of the distributions used to 



obtain them (see Figure 35). The dispersions of the distributions from which these Z and AY' are derived, a measure 
of the internal precision oi our photometry, are a ~ 0.01. 

ABSOLUTE CALIBRATION 

We placed our photometry in the Johnson- Cousins photometric system by observing 'Landolt| (|1992|) fields over a 
range of airmasses during the nights of our 2005 observing run. Since the first two nights were the most photometric 
of our run, we only used the standard fields observed those two nights to obtain the following transformations: 

B = 3.417 + 6 - 0.283X + 0.130(6 - v), (CI) 

V = 3.706 + V - 0.166X - 0.017(6 - v), (C2) 

/ = 3.065 + i - 0.049X - 0.020(v - i), (C3) 

where X is the airmass, 6, and i are the instrumental magnitudes in the instrumental system and 5, V, and / the 
corresponding magnitudes in the Johnson-Cousins photometric system. 

To fine-tune our calibration we decided to compare our photometry with Stetson] (|2000l) photometric standard stars. 
We found in some cases significant absolute offsets between Stetson's and our photometry. These zero points could 
be easily calculated (see Table 
Stetson's catalog, and only twef 



10 and Figure 36). However, only thirteen of our clusters had stars in common with 
ve had stars caHbrated in all three filters. 
To extend the comparison to more clusters, we looked for members of our sample already observed with at least 
one of the two HSTs optical wide-field imaging instruments, WFPC2 and ACS. We found that all of the clusters 



have already been observed at least in two of our filters, and most of them in all three (see Table 10). We retrieved 



the available data from the HST archive and de rived photometry from the stars i n these data u sm^ the programs 



HSTPHOT for the WFPC2 data ( |Dolphin||2000[ ) and DolPHOT for the ACS data (| Dolphin||2000[ ) . These programs 
have been specifically tailored to analyze data from these cameras. 

The comparison between the HST and ground-based data is complicated by the large resolution difference in the 
datasets and the highly crowded fields. When the offsets in t he p hotometries of individual stars of a cluster are plotted. 



there is a clear, asymmetric spread in the data (see Figure 37). Due to the higher resolution of the HST^ individual 
stars in the HST data are blends in the ground-based data, lb obtain the absolute offset between the two systems 
eliminating the blended stars effect we calculated a weighted mean of the offsets of the stars in common (see Figure 



37), where the weights were assigned according to the photometric errors and the dispersion in the offset distribution. 
The weighted mean was iteratively calculated, ignoring stars 2<j away from the mean value at each iteration. Whenever 
we have zero points calculated separately from independent WFPC2 and ACS data, a mean value was adopted. 
The results were then compared with the ones we found from the Stetson comparison, since our aim is to put all the 



photometry in the Stetson system (see Figure 38). A small offset for each filter was found and added to the individual 
cluster zero points {AB = -0.045 ± 0.012 mag, AV = -0.056 ± 0.007 mag, and AI = -0.025 ± 0.008 mag). These 
offsets may be due to the fact that the HST observations used in our absolute calibration are located closer to the 
center of the cluster than the Stetson stars used, or they may refiect small system zero point differences. The final 
offsets applied to the individual clusters can be seen in Table njl 
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Table 1 

Characteristics of the GCs in our sample according to the 2010 version of the [Ha?ri^ ( [19961 GGC 

catalog. 



Cluster name 


li 


b2 


Rgc^ 




[Fe 


/H] 


E(B 


-V) 


(m — 


M)v 


V 




NGC 6121 (M 4) 


350.97 


15.97 


5, 


.9 


2.2 


-1, 


.16 


0, 


,35 


12, 


.82 


-7, 


.19 


51.82 


NGC 6144 


351.93 


15.70 


2 


.7 


8.9 


-1, 


.76 


0, 


,36 


15, 


.86 


-6, 


.85 


33.35 


NGC 6218 (M 12) 


15.72 


26.31 


4, 


.5 


4.8 


-1, 


.37 


0, 


,19 


14, 


.01 


-7. 


.31 


17.28 


NGC 6235 


358.92 


13.52 


4, 


.2 


11.5 


-1, 


.28 


0, 


,31 


16, 


.26 


-6. 


.29 


11.18 


NGC 6254 (M 10) 


15.14 


23.08 


4, 


.6 


4.4 


-1, 


.56 


0, 


,28 


14, 


.08 


-7. 


.48 


18.47 


NGC 6266 (M 62) 


353.57 


7.32 


1, 


.7 


6.8 


-1, 


.18 


0, 


,47 


15, 


.63 


-9, 


.18 


11.28 


NGC 6273 (M 19) 


356.87 


9.38 


1, 


.7 


8.8 


-1, 


.74 


0, 


,38 


15, 


.90 


-9, 


.13 


14.57 


NGC 6287 


0.13 


11.02 


2 


.1 


9.4 


-2 


.10 


0, 


,60 


16, 


.72 


-7, 


.36 


6.96 


NGC 6304 


355.83 


5.38 


2 


.3 


5.9 


-0, 


.45 


0, 


,54 


15, 


.52 


-7, 


.30 


13.25 


NGC 6333 (M 9) 


5.54 


10.71 


1, 


.7 


7.9 


-1, 


.77 


0, 


,38 


15, 


.67 


-7, 


.95 


8.00 


NGC 6342 


4.90 


9.72 


1 
1 


7 


8.5 


-0, 


.55 


0, 


,46 


16, 


.08 


-D. 


.^z 


15.81 


NGC 6352 


341.42 


-7.17 


3 


.3 


5.6 


-0, 


.64 


0, 


,22 


14, 


.43 


-6. 


.47 


10.44 


iNvjV^ DOOO 






1, 


.4 


Q 9 


1 

-1 


.O ( 


n 
u, 


77 


1 7 


91 


-8, 


.07 


1 Ri 

lO.Ol 


NGC 6397 


338.17 


-11.96 


6, 


.0 


2.3 


-2 


.02 


0, 


,18 


12, 


.37 


-6, 


.64 


15.81 


NGC 6522 


1.02 


-3.93 


0, 


.6 


7.7 


-1, 


.34 


0, 


,48 


15, 


.92 


-7, 


.65 


15.81 


NGC 6541 


349.29 


-11.19 


2 


.1 


7.5 


-1, 


.81 


0, 


,14 


14, 


.82 


-8, 


.52 


13.04 


NGC 6553 


5.26 


-3.03 


2 


.2 


6.0 


-0, 


.18 


0, 


,63 


15, 


.83 


-7, 


.77 


7.66 


NGC 6558 


0.20 


-6.02 


1, 


.0 


7.4 


-1, 


.32 


0, 


,44 


15, 


.70 


-6, 


.44 


9.49 


NGC 6624 


2.79 


-7.91 


1, 


.2 


7.9 


-0, 


.44 


0, 


,28 


15, 


.36 


-7, 


.49 


18.97 


NGC 6626 (M 28) 


7.80 


-5.58 


2 


.7 


5.5 


-1, 


.32 


0, 


,40 


14, 


.95 


-8, 


.16 


11.23 


NGC 6637 (M 69) 


1.72 


-10.27 


1, 


.7 


8.8 


-0, 


.64 


0, 


,18 


15, 


.28 


-7, 


.64 


7.92 


NGC 6642 


9.81 


-6.44 


1, 


.7 


8.1 


-1, 


.26 


0, 


,40 


15, 


.79 


-6. 


.66 


9.77 


NGC 6656 (M 22) 


9.89 


-7.55 


4, 


.9 


3.2 


-1, 


.70 


0, 


,34 


13, 


.60 


-8, 


.50 


31.90 


NGC 6681 (M 70) 


2.85 


-12.51 


2 


.2 


9.0 


-1, 


.62 


0, 


,07 


14, 


.99 


-7, 


.12 


9.49 


NGC 6809 (M 55) 


8.79 


-23.27 


3 


.9 


5.4 


-1, 


.94 


0, 


,08 


13, 


.89 


-7, 


.57 


15.32 



^ Galactic longitude, in degrees. 

^ Galactic latitude, in degrees. 

^ Galactocentric distance, in kiloparsecs. 

^ Distance from the Sun, in kiloparsecs. 

^ Tidal radius, in arcminutes. 
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Table 2 

Summary of the observations with the Magellan telescope. 



Cluster 






UT Date 


Filter 


Exp. Time^ 


Airmass 


NGC6121 


16 23 40.0 


-26 30 29.8 


2005 Jun 2 


B 


3 X 1,3 X 25 


1 


37-1.42 










V 


3 X 1,3 X 20 


1 


43-1.48 










I 


3 X 1,3 X 15 


1 


50-1.57 


NGC6144 


16 27 18.7 


-26 00 28.4 


2005 May 31 


B 


3 X 5, 3 X 120 


1 


38-1.52 








V 


3 X 5, 3 X 120 


1 


55-1.65 










I 


1 X 1,3 X 3,3 X 75 


1 


77-1.92 


NGC6218 


16 47 18.6 


-01 55 51.6 


2005 Jun 1 


B 


3 X 5, 3 X 120 


1 


34-1.44 










V 


3 X 5, 3 X 120 


1 


46-1.55 










I 


1 X 1,3 X 3,3 X 75 


1 


56-1.66 


NGC6235 


16 53 29.8 


-22 09 37.8 


2005 May 31 


B 


3 X 5, 3 X 120 


1 


19-1.27 










V 


3 X 5, 3 X 120 


1 


28-1.34 










I 


1 X 1,3 X 3,3 X 75 


1 


36-1.44 


NGG6254 


16 57 13.0 


-04 04 58.3 


2005 Jun 1 


B 


3 X 5,3 X 120 


1 


14-1.16 










V 


3 X 5,3 X 120 


1 


17-1.20 










I 


1 X 1,3 X 3,3 X 75 


1 


21-1.25 


NGG6266 


17 01 17.2 


-30 05 42.8 


2005 May 30 


B 


3 X 5,3 X 120 


1 


01-1.02 








V 


3 X 5, 3 X 120 


1 


02-1.03 










I 


1 X 1,3 X 3,3 X 75 


1 


03-1.04 


NGG6273 


17 02 41.5 


-26 15 17.8 


2005 May 30 


B 


3 X 5,3 X 120 




1.00 








V 


3 X 5, 3 X 120 


1 


00-1.01 










I 


1 X 1,3 X 3,3 X 75 




1.01 


NGG6287 


17 05 13.8 


-22 41 28.2 


2005 Jun 2 


B 


3 X 5,3 X 120 


1 


24-1.32 










V 


3 X 5,3 X 120 


1 


33-1.42 










I 


1 X 1,3 X 3,3 X 75 


1 


44-1.53 


NGG6304 


17 14 37.1 


-29 26 43.2 


2005 May 31 


B 


3 X 5,3 X 120 


1 


06-1.08 








V 


3 X 5,3 X 120 


1 


09-1.11 










I 


1 X 1,3 X 3,3 X 75 


1 


12-1.15 


NGC6333 


17 19 16.1 


-18 30 00.0 


2005 Jun 2 


B 


3 X 5,3 X 120 


1 


02-1.03 










V 


3 X 5,3 X 120 




1.02 










I 


1 X 1,3 X 3,3 X 75 




1.02 


NGC6342 


17 21 14.5 


-19 34 13.1 


2005 May 31 


B 


3 X 5, 3 X 120 


1 


04-1.06 








V 


3 X 5, 3 X 120 


1 


02-1.04 










I 


1 X 1,3 X 3,3 X 75 




1.02 


NGC6355 


17 24 03.1 


-26 20 12.2 


2005 Jun 2 


B 


3 X 5, 3 X 120 


1 


03-1.06 










V 


3 X 5, 3 X 120 


1 


06-1.09 










I 


1 X 1,3 X 3,3 X 75 


1 


09-1.12 


NGC6352 


17 25 35.1 


-48 24 20.40 


2005 Jun 1 


B 


3 X 5, 3 X 120 




1.06 










V 


3 X 5, 3 X 120 


1 


06-1.07 










I 


1 X 1,3 X 3,3 X 75 




1.08 


NGC6397 


17 40 48.2 


-53 39 23.4 


2005 May 31 


B 


3 X 1,3 X 25 




1.11 








V 


3 X 1,3 X 20 




1.10 










I 


3 X 1,3 X 15 




1.10 


NGC6522 


18 03 38.5 


-30 00 59.6 


2005 May 30 


B 


3 X 5,3 X 120 


1 


01-1.02 








V 


3 X 5, 3 X 120 


1 


02-1.03 










I 


1 X 1,3 X 3,3 X 75 


1 


04-1.05 


NGC6541 


18 08 07.8 


-43 41 20.7 


2005 May 31 


B 


3 X 5,3 X 120 


1 


15-1.18 








V 


3 X 5, 3 X 120 


1 


18-1.22 










I 


1 X 1,3 X 3,3 X 75 


1 


23-1.26 


NGC6553 


18 09 20.1 


-25 53 26.6 


2005 May 31 


B 


3 X 5,3 X 120 


1 


02-1.03 








V 


3 X 5,3 X 120 


1 


03-1.06 










I 


1 X 1,3 X 3,3 X 75 


1 


07-1.09 


NGC6558 


18 10 23.1 


-31 44 48.0 


2005 May 30 


B 


3 X 5,3 X 120 


1 


07-1.10 








V 


3 X 5,3 X 120 


1 


11-1.13 










I 


1 X 1,3 X 3,3 X 75 


1 


14-1.18 


NGC6624 


18 23 45.2 


-30 20 39.7 


2005 May 30 


B 


3 X 5,3 X 100 


1 


21-1.24 








V 


3 X 5,3 X 100 


1 


25-1.30 










I 


1 X 1,3 X 3,3 X 60 


1 


31-1.36 


NGC6626 


18 24 37.3 


-24 51 10.8 


2005 Jun 2 


B 


3 X 5, 3 X 120 


1 


02-1.04 










V 


3 X 5, 3 X 120 


1 


04-1.05 










I 


1 X 1,3 X 3,3 X 75 


1 


06-1.07 


NGC6637 


18 31 28.0 


-32 19 55.6 


2005 Jun 2 


B 


3 X 5, 3 X 120 


1 


20-1.24 










V 


3 X 5, 3 X 120 


1 


25-1.30 










I 


1 X 1,3 X 3,3 X 75 


1 


17-1.19 


NGC6642 


18 31 58.7 


-23 27 36.0 


2005 Jun 2 


B 


3 X 5, 3 X 120 


1 


36-1.43 










V 


3 X 5,3 X 120 


1 


45-1.53 










I 


1 X 1,3 X 3,3 X 75 


1 


56-1.67 


NGC6656 


18 36 28.6 


-23 53 13.5 


2005 May 31 


B 


3 X 1,3 X 35 


1 


52-1.58 








V 


3 X 1,3 X 30 


1 


45-1.51 










I 


3 X 1,3 X 20 


1 


38-1.43 


NGG6681 


18 43 17.4 


-32 16 29.7 


2005 Jun 2 


B 


3 X 5,3 X 120 


1 


00-1.01 










V 


3 X 5,3 X 120 


1 


01-1.02 










I 


1 X 1,3 X 3,3 X 75 


1 


02-1.03 


NGG6809 


19 40 04.1 


-30 56 43.8 


2005 May 31 


B 


3 X 1,3 X 35 


1 


07-1.08 








V 


3 X 1,3 X 30 


1 


09-1.10 










I 


3 X 1,3 X 20 


1 


11-1.12 



^ Units of right ascension are hours, minutes, and seconds. 

Units of declination are degrees, arcminutes, and arcseconds. 
^ Number of exposures x exposure time, in seconds. 
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Table 3 

Summary of ACS observations from Snapshot program 10573. 



Cluster 




(^2000^ 


UT Date 


Filter 


Exp. Time^ 


NGC6218 


16 47 14.50 


-01 56 52.0 


2006 Feb 1 


^435 


120 










V555 


50 










^814 


20 


MGCd333 


17 19 11.80 


-18 30 59.0 


zOOd May 31 


^435 


2 X 340 










V555 


240 










^814 


90 


NGC6553 


18 09 17.50 


-25 54 28.0 


2006 Apr 4 


^435 


3 X 340 










V555 


300 










^814 


60 


NGC6624 


18 23 40.70 


-30 21 39.0 


2006 Jun 5 


^435 


360 










V555 


160 










^814 


65 


NGC6637 


18 31 23.20 


-32 20 53.0 


2006 Jun 6 


^435 


300 










V555 


120 










-^814 


50 



^ Units of right ascension are hours, minutes, and seconds. 

Units of declination are degrees, arcminutes, and arcseconds. 
^ Number of exposures x exposure time, in seconds. 



Table 4 

Summary of HST programs whose archival data 
have been used in our analysis. 



Cluster 


Camera 


Program 


Filters 


NGC6121 


WFPC2 


6116 


V555, ^814 




ACS 


10775 


Veoe, ^814 


NGC6144 


WFPC2 


11014 


^439, ^675 




ACS 


10775 


Veoe, ^814 


NGC6218 


WFPC2 


8118 


^439, V555 




ACS 


10775 


Veoe, ^814 


NGC6235 


WFPC2 


7470 


^439, ^555 




ACS 






NGC6254 


WFPC2 


6113 


Veoe, ^814 




ACS 


10775 


^606,-^814 


NGC6266 


WFPC2 


8118 


^439, V555 




WFPC2 


8709 


V555, -^814 




ACS 






NGC6273 


WFPC2 


7470 


^439,^555 




ACS 




NGC6287 


WFPC2 


6561 


^439, V555, /814 




ACS 






NGC6304 


WFPC2 


7470 


^439, V555 




ACS 


10775 


^606, ^814 


NGC6333 


WFPC2 


5366 


^555, ^814 




ACS 




NGC6342 


WFPC2 


7470 


^439, V555 




ACS 




NGC6352 


WFPC2 


5366 


V555, ^814 




ACS 


10775 


^606,-^814 


NGC6355 


WFPC2 


7470 


^439, V555 




ACS 






NGC6397 


WFPC2 


5929 


^439 , V555 




WFPC2 


7335 


V555, ^814 




ACS 


10775 


Veoe, ^814 


NGC6522 


WFPC2 


6095 


^439, ^555 




ACS 






NGC6541 


WFPC2 


5366 


^555, ^814 




ACS 


10775 


Veoe, ^814 


NGC6553 


WFPC2 


7307 


V555, ^814 




ACS 






NGC6558 


WFPC2 








ACS 


9799 


Veoe, ^814 


NGC6624 


WFPC2 


5324 


^439, V555 




WFPC2 


5366 


V555, ^814 




ACS 


10775 


Veoe, ^814 


NGC6626 


WFPC2 


6779 


V555, ^814 




ACS 






NGC6637 


WFPC2 


8118 


^439, ^55 




WFPC2 


5366 


V555, ^814 




ACS 


10775 


^606, ^814 


NGC6642 


WFPC2 


8118 


^439, ^555 




ACS 


9799 


Veoe, ^814 


NGC6656 


WFPC2 


5344 


^439, ^675 




WFPC2 


7615 


Veoe, ^814 




ACS 


10775 


^606,-^814 


NGC6681 


WFPC2 


8723 


^439, V555 




ACS 


10775 


Veoe, ^814 


NGC6809 


WFPC2 








ACS 


10775 


Veoe, -^814 
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Table 5 

Limits for the stars used in our 
analysis. 



Cluster 


VcL^ 






NGC 6121 


20.70 


1.9 




NGC 6144 


22.33 


1.7 




NGC 6218 


21.72 


1.8 




NGC 6235 


22.65 


2.1 


5.37 


NGC 6254 


22.07 


1.9 




NGC 6266 


21.51 


2.2 


6.84 


NGC 6273 


21.80 


2.3 


8.26 


NGC 6287 


21.90 


2.2 


4.13 


NGC 6304 


21.70 


2.9 


3.65 


NGC 6333 


22.29 


2.2 


5.59 


NGC 6342 


22.42 


2.3 


4.21 


NGC 6352 


21.86 


2.3 


5.61 


NGC 6355 


21.99 


2.2 


2.36 


NGC 6397 


20.80 


2.1 


8.18 


NGC 6522 


21.04 


3.4 


2.33 


NGC 6541 


22.30 


2.0 


8.63 


NGC 6553 


21.35 


2.4 


3.91 


NGC 6558 


21.48 


3.1 


1.81 


NGC 6624 


21.54 


2.7 


3.92 


NGC 6626 


21.04 


2.7 


4.95 


NGC 6637 


20.98 


2.2 


4.39 


NGC 6642 


19.99 


1.9 


2.22 


NGC 6656 


20.76 


2.4 




NGC 6681 


21.98 


2.2 


5.21 


NGC 6809 


21.89 


2.0 





^ V magnitude where the completeness 
limit is reached. 

^ Ay between the completeness limit 
and our dimmest observation for that 
cluster. 

^ Distance, in arcmin, from the center 
of the cluster where P{X = l|r) = 0.1. 
No distance given if the limit of the 
camera FOV is reached. 
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Table 6 

Paremeters used to obtain and fit the empirical |King| 
|l962l > models. 



Cluster 


r 












NGC 


6121 


1, 


,16 


51.82 


1069.8 


22.4 


14.2 


NGC 


6144 


0, 


,94 


33.35 


1067.8 


11.9 


36.5 


NGC 


6218 


0, 


,79 


17.28 


5336.6 


4.4 


8.1 


NGC 


6235 


0, 


,33 


11.18 


6140.2 


56.2 


82.3 


NGC 


6254 


0, 


,77 


18.47 


8905.6 


13.0 


13.3 


NGC 


6266 


0, 


,22 


11.28 


153639.3 


221.3 


222.8 


NGC 


6273 


0, 


,43 


14.57 


35673.5 


162.9 


154.8 


NGC 


6287 


0, 


,29 


6.96 


8549.7 


62.8 


66.1 


NGC 


6304 


0, 


,21 


13.25 


37443.9 


586.8 


548.4 


NGC 


6333 


0, 


,45 


8.00 


27214.7 


143.0 


125.8 


NGC 


6342 


0, 


,05 


15.81 


241666.3 


165.2 


145.6 


NGC 


6352 


0, 


,83 


10.44 


3042.5 


123.7 


177.5 


NGC 


6355 


0, 


,05 


15.81 


127282.0 


374.7 


650.0 


NGC 


6397 


0, 


,05 


15.81 


771500.0 


60.4 


40.5 


NGC 


6522 


0, 


,05 


15.81 


321899.2 


972.0 


1277.3 


NGC 


6541 


0, 


,18 


13.04 


190598.4 


85.3 


119.5 


NGC 


6553 


0, 


,53 


7.66 


11152.9 


428.5 


681.1 


NGC 


6558 


0, 


,03 


9.49 


376029.8 


609.4 


903.8 


NGC 


6624 


0, 


,06 


18.97 


232671.2 


310.2 


415.1 


NGC 


6626 


0, 


,24 


11.23 


57055.1 


377.1 


414.6 


NGC 


6637 


0, 


,33 


7.92 


11542.9 


115.4 


146.5 


NGC 


6642 


0. 


,10 


9.77 


4603.9 


111.9 


152.4 


NGC 


6656 


1, 


,33 


31.90 


4607.3 


159.8 


134.6 


NGC 


6681 


0, 


,03 


9.49 


1660360.2 


101.5 


118.6 


NGC 


6809 


1, 


,80 


15.32 


3253.0 


17.9 


17.0 



^ Core radius , in arcminutes, according to 2010 version of 
|Harris| |l996 | catalog. 

^ Tidal r adius, in arcminutes, according to 2010 version 
of Harris ( 1996 ) catalog. 

^ k constant from the fitted [King| ( | 196*21 profile. 

^ Density of non-member field stars trom the fit. 

^ Density of non-member field stars according to the Be- 

sangon model. 



Table 7 

Parameters used to obtain the Besangon model for the 
non-cluster stars in the observed field of NGC 6121, as an 
example of the technique. 



NGC 6121 



Field of view 

Field 



Extinction law 

Diffuse extinction 
Discrete clouds 

Selection on 

Intervals of magnitude 

Photometric errors 



Smah field 
1=350.97°; b=15.97° 
Solid angle=2.100 square degree 



An - 



0.0 mag/kpc 
=1.30; Distance^ 



:Opc 



15.04 <B < 24.48 

14.56 <V < 22.58 

13.18 < / < 22.53 
Error function: Exponential 
Band=B; A=0.006, B=21.75, C=0.861 
Band=y; A=0.006, B=22.24, C=0.901 
Band=/; A=0.007, B=23.66, C=0.997 
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Table 8 

Differential reddening, with respect to our ridgeline extinction 
zeropoints, of a selected set of coordinates across the field of our 
sampled clusters. Table fs] is published in its entirety in the 
electronic edition of A J. A~^portion is shown here for guidance 
regarding its form and content. 



^AE{B-V) Bandwidth^ 



NGC 6121 - M 4 



245.77881 


-26.63669 


-0.003 


0.009 


4.06 


245.77881 


-26.63571 


-0.004 


0.008 


4.03 


245.77882 


-26.63472 


-0.005 


0.008 


4.00 


245.77882 


-26.63373 


-0.007 


0.007 


3.97 


245.77883 


-26.63274 


-0.008 


0.006 


3.93 


245.77883 


-26.63175 


-0.010 


0.005 


3.90 


245.77883 


-26.63077 


-0.012 


0.002 


3.87 


245.77884 


-26.62978 


-0.012 


0.003 


3.84 



^ Units of right ascension are degrees. 
^ Units of declination are degrees. 

^ Units of the bandwidths used in the resolution map are arcmin- 
utes. 



Table 9 

Zero points to bring the photometry to the chip 2 
system. 



Chip Zb AY^ Zy AY{, Zi AYl 



1 
2 
3 
4 
5 
6 
7 
8 



-0.164 

0.0 
-0.107 
-0.767 
-0.010 
-0.137 
-0.268 
-0.212 



0.063 

0.0 
0.005 
0.048 
0.019 
0.010 
0.023 
0.055 



-0.160 

0.0 
-0.119 
-0.706 
0.013 
-0.173 
-0.260 
-0.169 



0.075 

0.0 
0.008 
0.043 
0.014 
0.006 
0.013 
0.042 



-0.036 

0.0 
-0.017 
-0.445 
0.083 
-0.118 
-0.185 
-0.028 



0.077 

0.0 
0.016 
0.078 
0.008 
0.021 
0.052 
0.044 
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Table 10 

Calibrating-star photometry available for the sampled clusters. 



Cluster Stetson^ WFPC22 ACS^ ACS-IOSTS^ if^'T missing^ Missing* 



NGC6121 


BVI 


VI 


VI 




B 


NGC6144 




B 


VI 






NGC6218 


BVI 


BV 


VI 


BVI 




NGC6235 




BV 






I 


NGC6254 


BVI 


VI 


VI 




B 


NGC6266 


BVI 


BVI 








NGC6273 


BVI 


BV 






I 


NGC6287 


BV 


BVI 








NGC6304 




BV 


VI 






NGC6333 




VI 




BVI 




NGC6342 




BV 






I 


NGC6355 




BV 






I 


NGC6352 


BVI 


VI 


VI 




B 


NGC6397 


BVI 


BVI 


VI 






NGC6522 


BVI 


BV 






I 


NGC6541 


BVI 


VI 


VI 




B 


NGC6553 




VI 




BVI 




NGC6558 






VI 




B 


NGC6624 




BVI 


VI 


BVI 




NGC6626 




BVI 








NGC6637 




BVI 


VI 


BVI 




NGC6642 




BV 


VI 






NGC6656 


BVI 


BVI 


VI 






NGC6681 


BVI 


BV 


VI 






NGC6809 


BVI 




VI 




B 



^ [stetson' ('2000') calibrating stars observed in these specific Johnson-Cousin filters. 

HST WFPC2 photometry from archival data available in the specified filters. The filters 
are /439'u;, f555w, fGOGw, and /814k;, and they are transformed by HSTPHOT to the 
Johnson-Cousins B, V, and /. 

^ HST ACS photometry from archival data available in the specified filters. The filters are 
f606w and /814k;, and they are transformed by DolPHOT to the Johnson-Cousins B, V, 
and /. 

^ HST ACS photometry from our project (SNAP 10573) available in the specified filters. 
The filters are /435k;, /555k;, and /814k;, and they are transformed by DolPHOT to the 
Johnson-Cousins B, V, and /. 

^ HST photometry for the cluster in the filter specified is missing. 

^ HST photometry and Stetson calibrating stars for the cluster in the filter specified are 
missing. 
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Table 11 

Offsets applied to transform to Johnson-Cousins system 
photometry. 



Offset in B Offset in V Offset in / 



Night 1 

NGC 6266 
NGC 6273 
NGC 6522 
NGC 6558 
NGC 6624 

Night 2 

NGC 6144 
NGC 6235 
NGC 6304 
NGC 6342 
NGC 6397 
NGC 6541 
NGC 6553 
NGC 6656 
NGC 6809 

Night 3 

NGC 6218 
NGC 6254 
NGC 6352 

Night 4 

NGC 6121 
NGC 6287 
NGC 6333 
NGC 6355 
NGC 6626 
NGC 6637 
NGC 6642 
NGC 6681 



-0.44 ± 0.02 

0.06 ± 0.02 
-0.07 ± 0.03 

0.01 ±0.02 



0.08 ± 0.01 
-0.02 ± 0.02 
-0.05 ± 0.02 



-0.02 ± 0.02 
-0.14 ± 0.03 
-0.03 ±0.03 
-0.20 ±0.02 
-0.13 ±0.02 
-0.27 ±0.03 
-0.64 ±0.03 
-0.13 ± 0.01 



-0.40 ± 0.02 
-0.10 ± 0.03 
-0.06 ± 0.01 
-0.06 ±0.01 
-0.03 ±0.01 



-0.13 ±0.02 
-0.05 ±0.02 
-0.02 ±0.01 
-0.13 ±0.03 
0.01 ± 0.01 
0.01 ± 0.02 
0.05 ±0.03 
0.02 ± 0.02 
0.04 ± 0.01 



-0.06 ± 0.01 
-0.02 ± 0.01 
-0.15 ± 0.02 



-0.38 ± 0.02 
-0.38 ± 0.03 
-0.05 ±0.02 
-0.09 ±0.02 
-0.21 ±0.03 
-0.46 ±0.02 
-1.25 ±0.02 
-0.25 ± 0.02 



-0.28 ± 0.02 
0.00 ± 0.03 
0.04 ± 0.02 
-0.05 ±0.01 
-0.04 ±0.01 



-0.10 ±0.02 
0.01 ±0.02 

-0.03 ± 0.01 
0.00 ± 0.02 
0.06 ±0.02 

-0.02 ± 0.02 
0.05 ± 0.02 



-0.01 ± 0.01 
-0.05 ± 0.01 
-0.15 ± 0.01 



-0.21 ± 0.02 
-0.05 ±0.02 
-0.17 ±0.02 

-0.46 ±0.02 
-0.20 ±0.02 
-1.42 ±0.02 
-0.10 ± 0.02 



-0.01 ±0.02 
-0.10 ±0.02 
0.02 ±0.02 
0.02 ±0.02 
0.01 ± 0.02 
0.05 ± 0.04 
0.16 ±0.07 
0.07 ± 0.02 
0.04 ± 0.02 



Note. — The numbers in itahc correspond to the straight 
comparison between the data from Stetson ( 2000 ) and our data. 
The numbers in standard type correspond to the comparison 
between our data and the HST data, plus the offset between 
HST data and Stetson's data. The errors reported are the 
dispersions of the distributions from which the offsets were 
calculated, which inform about the absolute precision of the 
photometry of individual stars. 
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Figure 1. Location, in galactic coordinates, of the clusters that follow the criteria to belong to our sample (solid triangles) and the rest 
of the inner GGCs (open triangles). The position of the 25 GGCs that we were able to observe and comprise our final sample is in red. 
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Figure 2. Radial density profile of the stars in the cluster from the observations (black) and from fitting a King profile (red), once 
the constant non-member field component has been subtracted. The straight green dashed line shows the non-member field component 
according to the fit. The straight blue dot-dashed line, showing the non-member field component according to the Besangon model, is 
plotted for comparison. The general disagreement in the inner regions between observations and models is due to the lower completeness 
of these regions, due to crowding effects. Whenever there are two black density profiles, one belongs to the space+ground observations, 
while the other corresponds to ground-only observations, and generally shows a lower completeness. 
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Figure 3. Probability to belong to the cluster as a function of distance to the center of the cluster. We should note, however, that both for 
the central regions where the crowding is more significant (as mentioned in Paper I), as for regions where extinction variations across the 
field are significant, the detection limit changes, in the case of the high crowding by being lower, and in the case of significative extinction 
by, although being nominally the same, corresponding to different intrinsic luminosities in different parts of the field. This is going to 
produce deviations from the model as we have seen in Figure f2| In our analysis, we assume that we are missing an equal percentage of 
stars from the cluster and from the field due to these facts, anonence we use the values provided in this figure. 
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Figure 4. Observed dereddened CMD for M4 at four different distances from the cluster center, with the probabihty of every star to 
belong to the cluster as a function of position in the sky, color and magnitude, P(X = l|r, c, m), represented by the different colors of the 
stars as indicated in the color bar. Notice that we restrict our analysis to stars brighter than the completeness limit shown in Table [s] 
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Figure 7. On top, CMDs of the cluster NGC 6121 - M 4, before and after being differentially dereddened. Only our Magellan photometry 
was used to build the B — V vs. V CMD. ACS photometry (from project 10775) and Magellan photometry were used to build the V — I vs. 
V CMD. Color bars show the range in the densities of stars in the CMD (xlO^ stars per square magnitude). On the bottom, extinction 
map (left) for the cluster field, along with its resolution (middle) and its precision (right), as provided by our technique. The x marks the 
position of the cluster center. The color code gives, respectively, the color excesses E(B — V) for the extinction map, the bandwidths used 
in the resolution map, and the standard deviation cr of the color excesses in the precision maps. 
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Figure 8. As in Figure [T] but for the cluster NGC 6144. Only our Magellan photometry was used to build the B - V vs. V CMD. ACS 
photometry (from project 10775) and Magellan photometry were used to build the V — I vs. V CMD. Notice that the B — V vs. V CMD 
could not be correctly calibrated in color using the method described in the text because of the lack of calibrating data in the B filter. 
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Figure 9. As in Figure [7| but for the cluster NGC 6218 - M 12. ACS photometry (from our project 10573) and Magellan photometry 
were used to build the B — V vs. V CMD. ACS photometry (from project 10775) and Magellan photometry were used to build the V — I 
vs. V CMD. 
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Figure 10. As in Figure ItI but for the cluster NGC 6235. Only our Magellan photometry was used to build the CMDs in both colors. 
Notice that the V — I vs. vCMD could not be correctly calibrated in color using the method described in the text because of the lack of 
calibrating data in the / filter. 
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Figure 11. As in Figure|7| but for the NGC 6254 - M 10. Only our Magellan photometry was used to build the B-V vs. V CMD. ACS 
photometry (from project 10775) and Magellan photometry were used to build the V — I vs. V CMD. 
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Figure 12. As in Figure ItI but for the cluster NGC 6266 - M 62. Only our Magellan photometry was used to build the B — V vs. V 
CMD. WFPC2 photometry (from project 8709) and Magehan photometry were used to build the V - I vs. V CMD. 
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Figure 13. As in Figure [T] but for the cluster NGC 6273 - M 19. Only our Magellan photometry was used to build the CMDs in both 
colors. 
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Figure 14. As in Figure [7| but for the cluster NGC 6287. Only our Magellan photometry was used to build the CMDs in both colors. 
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Figure 15. As in Figure [T] but for the cluster NGC 6304. Only our Magellan photometry was used to build the S - V vs. V CMD. ACS 
photometry (from project 10775) and Magellan photometry were used to build the V — / vs. V CMD. 
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Figure 17. As in Figure ItI but for the cluster NGC 6342. Only our Magellan photometry was used to build the CMDs in both colors. 
Notice that the V — I vs. vCMD could not be correctly calibrated in color using the method described in the text because of the lack of 
calibrating data in the / filter. 
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Figure 18. As in Figure [T] but for the cluster NGC 6352. Only our Magellan photometry was used to build the B-V vs. V CMD. ACS 
photometry (from project 10775) and Magellan photometry were used to build the V — I vs. V CMD. 
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Figure 19. As in Figure f7| but for the cluster NGC 6355. Only our Magellan photometry was used to build the CMDs in both colors. 
Notice that the V — I vs. vCMD could not be correctly calibrated in color using the method described in the text because of the lack of 
calibrating data in the / filter. 
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Figure 20. As in FigurejTl but for the cluster NGC 6397, before and after being differentially dereddened. Only our Magellan photometry 
was used to build the B -^V vs. V CMD. ACS photometry (from project 10775) and Magellan photometry were used to build the V — I 
vs. V CMD. 
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Figure 21. As in Figure [7| but for the cluster NGC 6522. Only our Magellan photometry was used to build the CMDs in both colors. 




Figure 22. As in Figure[7| but for the cluster NGC 6541. Only our Magellan photometry was used to build the B -V vs. V CMD. ACS 
photometry (from project 10775) and Magellan photometry were used to build the V — I vs. V CMD. 



47 





I I I I I I I I I I I I I I I i I I I I i I I I I I I I I ■-0.1 I I I I I I I I I I I I I I I i I I I I i I I I I I I I I ■ I I I I I I I I I I I I I I I i I I I I i I I I I I I 

272.2 27225 2723 27235 2724 27245 2722 27225 2723 27235 2724 27245 2722 27225 2723 27235 2724 27245 

Right Ascension Right Ascension Right Ascension 



Figure 23. As in Figure m but for the cluster NGC 6553. ACS photometry (from our project 10573) and our Magellan photometry were 
used to build the CMDs inboth colors. 
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Figure 24. As in Figure [tI but for the cluster NGC 6558. Only our Magellan photometry was used to build the B -V vs. V CMD. ACS 
photometry (from project 9799) and Magellan photometry were used to build the V — I vs. V CMD. Notice that the B — V vs. V CMD 
could not be correctly calibrated in color using the method described in the text because of the lack of calibrating data in the B filter. 
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Figure 25. As in Figure [t] but for the cluster NGC 6624. ACS photometry (from our project 10573) and Magellan photometry were 
used to build the B — V vs. V CMD. ACS photometry (from project 10775) and Magellan photometry were used to build the V — I vs. 
V CMD. 
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Figure 26. As in Figure [T] but for the cluster NGC 6626 - M 28. Only our Magellan photometry was used to build the B — V vs. V 
CMD. WFPC2 photometry(from project 6779) and Magehan photometry were used to build the V - 1 vs. V CMD. Notice that the B-V 
vs. V CMD could not be correctly calibrated in color using the method described in the text because of the lack of calibrating data in the 
B filter. 
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Figure 27. As in Figure m but for the cluster NGC 6637 - M 69. ACS photometry (from our project 10573) and Magehan photometry 
were used to build the B vs. V CMD. ACS photometry (from project 10775) and Magellan photometry were used to build the V — I 
vs. V CMD. 
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Figure 28. As in Figure [T] but for the cluster NGC 6642. Only our Magellan photometry was used to build the B-V vs. V CMD. ACS 
photometry (from project 10775) and Magellan photometry were used to build the V — I vs. V CMD. 
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Figure 29. As in Figure |7| but for the cluster NGC 6656 - M 22. Only our Magellan photometry was used to build the B — V vs. V 
CMD. ACS photometry (from project 10775) and Magellan photometry were used to build the V — I vs. V CMD. 
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Figure 30. As in Figure |7| but for the cluster NGC 6681 - M 70. Only our Magellan photometry was used to build the B — V vs. V 
CMD. ACS photometry (from project 10775) and Magellan photometry were used to build the V — I vs. V CMD. 
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Figure 31. As in Figure |7| but for the cluster NGC 6809 - M 55. Only our Magellan photometry was used to build the B — V vs. V 
CMD. ACS photometry (from project 10775) and Magellan photometry were used to build the V — I vs. V CMD. 
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Figure 32. Comparison of the E(B-V) extinction maps provided by our technique (left) with the SFD maps (right). We have lowered the 
resolution of our maps (see text) to compare them more easily with the SFD maps. From the comparison of both maps for every GC we 
have obtained a reddening zero point for our map, that we have added to facilitate the comparison. 
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Figure 33. Comparis on of the E(B-V) extinction maps provid ed by our technique (left) with t he extinction maps (right) pro vided by 
[Mochejskaet al.'f2002^ for NGC 6121, by von Braun et al."f2002^ for NGC 6218 and NGC 6254, by'Gerashchenko & Kadla f^OOT) for NGC 
d266, and by Meitsch ^ Richtler ( 1999) for INGU 6342 and NGU 6553. Ah the maps from the hterature show relative ditterential extinctions 
with respect to the extinction zero point of a fiducial region, except for NGC 6342 and NGC 6553, that show absolute extinctions. Notice 
that the extinction zero points between our maps and those from the literature can differ. 
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Figure 34. Pointings used to make the instrumental calibration, bringing the photometry from the different chips to the chip2 system. 
We chose these pointings because they allowed us to observe stars in common between chip2 and all the other chips in just 3 pointings 
(chips 1, 5 and 6 have stars in common with chip 2 in the left configuration and chips 3, 4, 7, and 8 have stars in common with chip2 
in the right configuration). We used stars from the cluster NGC6397, one of the clusters in our sample, to do this calibration. The blue 
pointing was the same as the one originally used in our observing run. That way we were able to compare the photometries of our original 
observing run and the calibrating run for all the chips in the camera and take care of any effect produced by a change in the sensitivity of 
the chips (see text). 
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Figure 35. Calculation of the zero points to bring the photometry to the chip2 system. In this figure we have plotted, as an example of 
the technique, the different zero points that we have to obtain to move the V photometry from the chip3 to the chip2 reference system (see 
equations in the text). The solid lines show the final average values adopted. 
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Figure 36. Comparison of our raw ground-based photometry with the |Stetson| ([2000') cahbrating stars for the clusters and filters available. 
We plot the magnitudes of the stars in our photometry for the different hlters versus Stetson's values minus ours. The lines show where 
the average offsets lie. 
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Figure 37. Comparison of our raw ground-based photometry with the HST photometry. In this figure we show the steps to calculate 
the offsets in NGC 6218, as an example of the technique. We plot the magnitudes obtained from our ground-based observations versus the 
HST photometry minus our ground-based photometry. We can observe the existence of the spread mentioned and explained in the text. 
The average offsets, shown by the lines, are calculated using a clipped weighted average (see text). 
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Figure 38. Comparison of our photometry corrected with the available ff5'T offsets, with the |Stetson| ( |2QQQ] ) calibrating stars. We plot 
the average differences of Stetson's values minus ours, for every cluster and filter. The lines show the weighted average offsets for every 
filter. 



